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ABSTRACT 

We investigate the properties of massive galaxies at z 1-3.5 using Hubble Space Telescope observations 
at optical wavelengths, ground-based near-infrared (IR) imaging, and Spitzer Space Telescope observations at 
3-24 /im. From /Tj-selected galaxies over a ~ 130 arcmin^ field in the southern Great Observatories Origins 
Deep Surveys (GOODS-S), we identify 153 distant red galaxies (DRGs) with (7-/rs)vega > 2.3. This sample 
is approximately complete in stellar mass for passively evolving galaxies above 10" M.q and z < 3. Of the 
galaxies identified by this selection, roughly half are objects whose optical and near-IR rest-frame light is 
dominated by evolved stars combined with ongoing star formation (at ZmeA ~ 2.5), and the others are galaxies 
whose light is dominated by heavily reddened (A1600 ^4-6 mag) starbursts (at Zmed 1-7). Very few of the 
galaxies (< 10%) have no indication of current star formation. The total star-formation rates (SFRs) including 
the reradiated IR emission for the DRGs are up to two orders of magnitude higher than those derived from the 
UV luminosity corrected for dust reddening. We use population synthesis models to estimate stellar masses 
and to study the stars that dominate the rest-frame UV through near-IR light in these galaxies. DRGs at 
z ~ 1 .5 - 3 with stellar masses > 10" M.q have specific SFRs (SFRs per unit stellar mass) ranging from 0.2 to 
10 Gyr~', with a mean value of ^ 2.4 Gyr"'. Based on the X-ray luminosities and rest-frame near-IR colors, 
as many as one-quarter of the DRGs may contain AGN, implying that the growth of supermassive black holes 
coincides with the formation of massive galaxies at z > 1.5. The DRGs with > lO" M.q at 1.5 < z < 3 
have integrated specific SFRs greater the global value over all galaxies at this epoch. In contrast, we find that 
galaxies at z 0.3-0.75 with M > 10" M.q have integrated specific SFRs less than the global value, and 
more than an order of magnitude lower than that for massive DRGs atz~1.5-3. Atz<l, lower-mass galaxies 
dominate the overall cosmic mass assembly. This suggests that the bulk of star formation in massive galaxies 
occurs at early cosmic epochs and is largely complete by z 1.5. Further mass assembly in these galaxies 
takes place with low specific SFRs. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: high- 
redshift — galaxies: stellar content — infrared: galaxies 
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L INTRODUCTION 

Most of the stellar mass in galaxies today apparently 
formed during the relat ively short period between z ~ 
3 and 1 (e .g. Dickinson et alJ l2003t jRudnick et al. 200^ 
'Fontana et al. 2004; Glazebrook et al.''2004). Some early- 
type galaxies appear as soon as z ^ 1.5-2 ( Dunlop 
1996: jSt?inrad et alJIlQQTtlMcCarthv et a l. 2004; Daddi 
12005 alK By z ~ 1, there is a significant population 
of galaxies with red colors and morphologies consis- 
tent with passively-evolving early-type galaxies, imply- 



ing t hey formed their stel l ar population s at Zform 



> 2-3 

(e.g., IC imatti et al. 2002a; Moustakas & Somerville 200l 
^vloustakas et al. 2004; Papovich et al. 2005 ; Treu et al,.200a 
McCarthv 2004^ for a review). The cosmic star-formation 
rate (SFR) density has declined by roughly a factor 10 from 
z ~ 1 to the present-day (e.g. Hopkins 2004, and refer- 
ences therein). During the time since z 1 the stellar mass 
in passively-evolving, early-type galaxies has increased by 
less th an a factor 2 (e.g.|Brinchmann & Ellis 2000; Bell et aj] 
120041) , and at the present epoch, roughly one-third of all stars 
exist in such galaxies (Baldrv et al. 2004). 

Although massive galaxies appear to have formed most 
of their stars at epochs prior to z 1 - 2, we have few 
constraints on how they assembled their stellar mass. One 
hypothesis is that galaxies "downsize" and star formation 



2 



PAPOVICH ET AL. 



shifts from galaxies at the high to the lo w end of the mass 
function with decre a sing redshifts (e.g., Cowie et all 119991: 
Fontana et all |2003| iHeavens et al.1 120041: iKauff mann et all 



2004; Bau er et all l2005t ICaputi et all l2005n juneau et all 

2£)05; Perez-Gonzalez et al. 2005). Another possibility is that 
massive galaxies assemble their stellar mass early-on, either 
in situ, "closed-box " formation events with subsequent pas - 
sive evolution (e.g.. 'Egg en. Lvnden-Bell. & Sandag^l 19621) . 
or in low-mass systems, which then coalesce to form 
large galaxies with little subsequent star formatio n (e.g., 
paugh^etal. .l998; Kauffm ann & Chariot 1998: Cimatti et alJ 

Testing these proposals has been frustrated by difficulties 
in conducting a complete census of star-forming galaxies 
and massive galaxies at z > 2, when such systems are ex- 
pect ed to experience a peak in their stellar assembly rates 
(e.g. iDe Lucia et al.l 2005: Nagamine et al. 2005). Ultravi- 
olet (UV)-luminous star-forming galaxies at these redshifts 
are identified by the characteristic "break" in their colors due 
to neutral hydrogen absorption that attenuates the flux short- 
ward of Lyman a (1216 A) and the 912 A Lyman limit (e.g. 
ISteidel et allT99 6: Giavalisco 2002, for a review). These UV- 
dropout, Lyman-break galaxies (LBGs) dominate the UV lu- 
minosity density at z 2- 6, and possibly the global SFR den- 
sity at these redshifts (e.g.. Steidel et al 1999: Bouwe ns et"al] 
12004: Giavalisco et al. 2004b). However, the UV-dropout 
technique is primarily sensitive to galaxies with ongoing, rela- 
tively unreddened star formation. Surveys with SCUBA have 
identified a population of sub-mm galaxies at z > 2 that emit 
the bulk of their bolometric luminosity at infrared (IR) wave- 
lengths (see Blai n et al. 2002. for a review), and may con- 
tribute substantially to the cosmic SFR (Bargeret al. 2000). 
Their inferred space densities and SFRs suggest that they 
could be the progenitors of the most mass ive present-day 
galaxies (Blain et al. 2004; Chapman et al."2005). Neither the 
LBG nor the sub-mm-galaxy populations necessarily provide 
a full sample selected by stellar mass, so how they participate 
in the formation of present-day massive galaxies is poorly 
understood. Hence, it is unclear what fraction of present-day 
galaxies pass through such stages during their assembly. 

Theoretically, the star-formation histories of massive 
galaxies are also poorly known. Galaxies within mas- 
sive halos have short cooling times and tend to convert 
all their gas into stars rapidly, unless feedback is in- 
voked from stellar winds and supernovae to reheat the 
gas o r prevent it from cooling (e.g. Hernquist & Spring^ 
I2;003t Springel et al. 2005a b). Models predict that mas- 
sive galaxy halos continue to accrete smaller aggregates 
up to the current epoch, which rejuvenates star-formation 
and predicts galaxy colors that are too blue compared 
to observations (e.g. Somerville, Primack, & Faber 200^ 
unless very large dust extinctions or non-sta ndard stel- 
lar in itial mass functions (IMFs) are invoked fcBau^h et alJ 
l20()5t Nagami ne et alJ |2005). Some recent theoretical 
models suppress star-formation at late times in massive 

Eies using feedback from active-galactic nuclei (A GN; 
itoetal. 2001; Di Matteo et al. 2005; Hopki ns e^tal] 
Springel et al., ,2005a) . This process may pro- 
vide the impet us for the present-day black-hole-bulge- 
mass correlation ('Magorr ian et alJl998H Gebhardt et alJ200(t 
[Kauffmann et al. 2004). We require observations of star- 
formation and AGN activity in high-redshift, massive galax- 
ies to improve our understanding of how such systems form. 



Surveys using deep, near-IR observations have iden- 
tified high-redshift, n iassive-galaxy candida t es with red 
near-IR colors (e.g., 'Dickinson et al.l 120001 iT otani et alJ 
12001; Imetal. 2002; Daddi et al. 20041 lYan et al. 2004a). 
TVanxetal. (2003) used deep JHKs observations from 
VLT/ISAAC to identify a population of galaxies with 
(J-Ks)ysga > 2.3 mag in the Faint IR Extragalactic Survey 
(FIRES). In principle, this color selection identifies galaxies 
that have a strong Balmer/4000 A break between the J and 
Ks bands at z ~ 2-3.5, down to an approximately complete 
limit in stellar mass. Subsequent analysis has concluded that 
these distant red galaxies (DRGs) are m ostly massive, old, 
and actively forming stars at z ~ 1.5-3.5 ('van Dokkun i et all 
2003; Forster-Schreiber et a l] 120041: Rubin et al. t 2004t 
Knud sen et alJ 120051: IReddv et alJ I2005ID although some 
appear to be co mpletely devoid o f star-formation and 
passively evolving (Lab be et alJl2005 |). while other s appear 
to host powerful AGN (Ivan Dokkum et alJ 12004ft . The 
inferred stellar masses of DRGs at z 2 - 3 are similar to 
those of local early-type galaxies (Forster-Schreibe r et alJ 
120041) . They are generally higher th an those inferred from 
LBG samples at similar redshifts dSawicki & Ye 
iPapovich. Dickinson. & Ferguson*2001l lShaplev et alJ2001h . 
although some overlap between the two clearly exists 
iShaplev et al. 2005; Reddy et al. 2005). The estimated stel- 
lar population ages of the DRGs suggests that they have bee n 
forming stars since z ~ 5 - 6 ( Forster-Schre iber et a 
Thus DRGs may represent the older stellar populations 
formed in higher-redshift LBGs (e.g. Papovich et al. 2004a), 
and these galaxies possibly lin k the UV -luminous LBGs to 
the sub-mm galaxies ( van Dokkum et al. 2004). 

In addition, extinction may con tribute to the red colors 
of some of the galaxies. ISmail et aL (k2002) find that dust- 
extincted starbursts at z ~ 1-2 selected as extremely red ob- 
jects (EROs) typically have red J-K colors similar to DRGs. 
Using current hierarchical model predictions, Nagamine et aD 
(2005) suggest that massive galaxies in formation at z > 1 
may be heavily reddened by dust. If so, then most of their 
emission should appear in the thermal IR. In fact, up to half 
of EROs (selected with red R-K or /?-[3.6/im] colors) 
at z > 1 are detected in the thermal IR by Spitzer/MlPS at 
24 fim (Wilson etal. 2004; Yan et al. 2004b). In addition, 
the inferred evolution of the galaxy population responsible 
for the IR number counts implies a substanti al population 
of galaxies at z ~ 1-3 (Papovich et al. 2004b: Caput i et all 
^05; Perez-Gonzalez et al. 2005), and the evolution in the 
IR luminosity function suggest that luminous IR galaxies 
(LIRGs, LiR = 10"~'^ L0), and ultraluminous IR galaxies 

i LiR > 10'^ L(n) dominate the IR luminosity densi ty at z > 1 
Le Floc'h e t al.'2005l: IPerez-Gonzalez et alJ2005h . 

For this study, we selected galaxies at z ~ 1-3.5 with 
(7-/:s)vega > 2.3 mag {[J-K,]ab > 1-37 mag) from a /ri- 
band selected catalog in the southern Great Observatories 
Origins Deep Survey (GOODS-S) field. At the magnitude 
limit of Ks < 23.2 [AB], this sample is approximately com- 
plete in stellar mass for passively evolving galaxies with M > 
10" A^o for z < 3. We also use IR and X-ray observations 
to constrain the star-formation and AGN processes in these 
galaxies. The AGN connection in these galaxies is explored in 
more detail in a forthcoming paper (L. A. Moustakas et al., in 
preparation). In §|2l we summarize the data and we define the 
galaxy samples. In § |3l we describe broad properties of the 
DRGs using their rest-frame UV to near-IR colors. In § |3 
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we compare the IR luminosities for the galaxies with those 
derived from their UV luminosity and measured extinction, 
and we compare the DRGs to other galaxies at 1 .5 < z < 3.5. 
In § |5l we use population synthesis models to estimate the 
properties of the galaxies' stellar populations. In §|6j we dis- 
cuss the relationship between the SFRs and stellar population 
properties for the ensemble of galaxies, and we comment on 
the presence of AGN. We also compare the distribution of 
SFR as a function of galaxy stellar mass at high redshift, and 
we compare with results at lower redshifts (z ^ 0.3-0.75). In 
§0we present our conclusions. 

Throughout this paper we use a cosmology with Oxotai = 1, 
r^M = 0.3, A = 0.7, and Hq = 70 km s"' hjo Mpc"' where 
hjo = 1. Unless otherwise specified, we present all magni- 
tudes in the AB system, hjab = 23.9-2.51og(/y/l fiiy). We 
denote galaxy magnitudes from the HST ACS bandpasses 
F435W, F606W, F775W, and F850LP as B435, Veoe, his, 
and Z850, respectively. Similarly, where applicable we de- 
note magnitudes from the HST WFPC2 and NICMOS band- 
passes F300W, F450W, F606W, F814W, Fl lOW, and F160W 
as t/300, -B450, V6O6, hi4, -/no, and Hia), respectively. We also 
denote magnitudes from the four Spitzer IRAC channels as 
[3.6/im], [4.5/xm], [5.8/im], and [8.0/im], respectively. 

2. THE DATA AND SAMPLE DEFINITIONS 

The GOODS-S field center is located in the southern Chan- 
dra Deep X-ray field (CDF-S) at 3'>32™30\ -27°48'20", 
which provides 1 Ms of imaging in the soft (0.5-2 keV) and 
hard (2-8 keV) X-ray bands in this field. The available 
observations include imaging with HST/ACS in four pass- 
bands, B435, V6O6, '775, and zg5o over 160 arcmin^, ground- 
based near-IR imaging from VLT/ISAAC in the JK, bands 
over 130 arcmin^ (with //-band imaging over ~ 50 arcmin^), 
and IR imaging from SpitzerflRAC in four bands, [3.6/im], 
[4.5/im], [5.8yum], [8.0Aim]. 

The ACS obse rvations and data reduction are described in 
iGiavali sco et al. (2004a). The images have a PSF FWHM 
~ 0."125, and the limiting lOcr sensitivities are B435 = 27.8, 
Veoe = 27.8, ijjs = 27.1, zggo = 26.6, measured in 0"2-diameter 
circular apertures. The Spitzer/IRAC images have a PSF 
FWHM ranging from ~ 1 "5 at 3.6 pm to ~ 2"at 8 pm, and for 
isolated point sources achieve Sa limiting sensitivities rang- 
ing from 0.11 fi}y at 3.6 pm to 1.66 p]y at 8 pm (M. Dickin- 
son, et al., in preparation). 

The near-IR ISAAC imaging is from the version 1.0 re- 
lease (B. Vandame et al., in preparation), available on the 
ESO/GOODS webpages."* The ISAAC data have excellent 
image quality (full-width at half maximum, FWHM w 0."45) 
with mean exposures times of 14000 and 24000 s in 7 and 
Ks, respectively, reaching limiting magnitudes of 7 = 24.7, 
H = 24.1 and K, = 24.1 (10 a) in l"-diameter apertures, al- 
though the depth varies over the GOODS-S field. 

Spitzer imaged the CDF-S field with MIPS at 24, 70, 
and 160 pm under 5/?/rzer/Guaranteed Time Observer (GTO) 
time. Here, we focus exclusively on the 24 pm imaging, 
which was reduced using the instrument team Data Analysis 
Tool ( Gordon et al. 2005). The GTO MIPS imaging covers 
1° X 0.5° with a FWHM « 6 arcsec, and covers all of the 
GOODS-S field and most of the ESO imaging and COMBO- 
17 surveys (see § I2.5> . 

2.1. GOODS Source Cataloging and DRG Sample Selection 
http://www.eso.org/science/goods/releases/20040430/ 



We use a source catalog selected from the ISAAC /Tj-band 
data. We rebinned the ACS data to the pixel scale of ISAAC, 
and convolved the ACS images to match the image quality of 
the ISAAC images. Source catal ogs were the n constructed us- 
ing the SExtractor software (Berlin & Arnouts 1996) by first 
locating sources on the /Ts-band image, then measuring pho- 
tometry in matched apertures on the HST/ACS and ISAAC 
images. We measured photometry in each band in isophotal 
apertures defined from the /Tj-band image. We then scaled 
these to total magnitudes using the difference between the 
/Tj-band isophotal magnitude (SExtractor MAG_ISO) and the 
/Tj-band magnitude measured in a "total", elliptical aperture 
defined by the Kron radius (SExtractor MAG_AUTO). Pho- 
tometric uncertainties are derived by SExtractor after adjust- 
ing the image rms maps to account for the correlated noise 
properties introduced by drizzling.'^ The SExtractor-derived 
uncertainties still likely underestimate the true errors, because 
they do not account for systematic errors in the measurements 
themselves. Therefore, we have included an additional error 
of cTsys/Ziy ~ 3%, added in quadrature to the uncertainties on 
the ACS and ISAAC photometry. 

We detected objects in the SpitzerflRAC images using a 
weighted-sum image of IRAC channels 1 and 2. Magnitudes 
were then measured in each IRAC band in 4"-diameter aper- 
tures, and we applied aperture corrections of 0.30, 0.34, 0.53, 
and 0.67 mag to the bands [3.6/im], [4.5/im], [5.8/im], and 
[8.0/im], respectively, to correct to total magnitudes. We did 
not attempt to measure photometry on versions of the ISAAC 
data PSF-matched to the IRAC data, because of the signifi- 
cantly poorer image quality of the IRAC data. The aperture 
corrections are based on Monte Carlo simulations using artifi- 
cial, compact sources added to the real images, and are appro- 
priate for sources with half-light radii < 0"5, such as those of 
interest here. These simulations also allow us to estimate the 
error on the IRAC photometry as a function of flux density. 
The uncertainties are similar for IRAC channels 1 and 2, and 
range from Am « 0.03 mag at m w 21 mag to Am « 0.3 mag 
at m « 25 mag. For IRAC channels 3 and 4, they range 
from Am « 0.05 mag at m « 21 mag to Am « 0.4 mag at 
m « 25 mag. 

We used {J-Ks)^^,^^ > 2.3 mag to identify DRGs. In the 
GOODS-S field, we found 153 of them to a signal-to-noise 
ratio S/N(Ks) > 10 limit, which is the median S/N of ob- 
jects at Ks < 23.2 mag within the "total", MAG_AUTO aper- 
tures, and with J-K, colors measured in the seeing-matched 
MAG_AUTO apertures. Restricting sources to the S/N re- 
quirement is appropriate as the ISAAC depth varies over 
the GOODS-S field. Using a S/N limit also ensures that 
we can derive robust colors from the ACS and ISAAC data, 
greatly improving the accuracy of our SED modeling and pho- 
tometric redshifts. We inspected the ACS image for each 
DRG at their native resolution to identify objects resulting 
from chance galaxy-galaxy alignments along the line of sight. 
In one case the DRG does appear to involve multiple ACS 
sources, blended at the K-band resolution. We exclude this 
object although its inclusion does not effect the results in this 
paper A unique IRAC source exists for 132 of the 153 DRGs 
with a matching radius of r < 0"5. We visually inspected each 
object to verify that the matched IRAC source corresponds to 
the DRGs in the /Tj-band image. The unmatched objects suf- 
fer from crowding in the IRAC images from other sources 
within ~ 1-2". The IRAC flux from these non-detected 

see |ttp://arcMve.stsci.edWpub/Msp/goods/vl/h_goo3i_vn()~rdni.htnil] 
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Fig. 1. — Distribution of redshifts for tlie DRG and HDF-N galaxy sam- 
ples. The top panel shows the redshift distribution of DRGs with no 24 /im 
detections. The middle panel shows the redshift distribution of DRGs with 
/y(24^m) > 50 /iJy. The bottom panel shows the redshift distribution of the 
24 galaxies from the HDF-N with 1.5 < z < 3.5 and K, < 23.2 mag. The 
dashed lines indicate the redshift limits of the sample. The shaded histogram 
shows the distribution for those HDF-N galaxies with /i,(24/^m) > 10 fj,]y. 



DRGs is either completely blended within the isophote of the 
neighbor, or confusion with the neighbor offsets the centroid 
of the IRAC flux past the matching criterion. 

Spectroscopic redshifts are available for 12 of the DRGs 
(Szokoly et al. 2004; Mignoli et al. 2005; Vanzefla et al. 
2005; D. Stern 2005, private communication). We there- 
fore supplement the redshift information using photometric 
redshifts from Mobasher et al. (2004 and in preparation). 
The high-quality dataset allows for accurate photometric red- 
shifts. The median photometric-redshift uncertainty for all 
^Tj-band sources with spectroscopic re dshifts is fe /(l +z) ~ 
0.1, similar to the accuracy found by C aputi et al. ' (2005) us- 
ing similar data in the same field. For the 12 DRGs with 
spectroscopic redshifts, the photometric-redshift accuracy is 
even better, ~ 0.04. A photometric redshift is 

unavailable for one of the DRGs, whose ACS and ISAAC 
7-band photometry have S/N< 1. We ignore this DRG for 
any analysis requiring a redshift (e.g., luminosities, stellar 
masses, etc.). Figure ^ shows the redshift distribution. The 
redshifts range from z — 0.8-3.7, with a median redshift 
(z) = 2.2. This is in broad agreement with the re dshift distri- 
bution for DRGs reported in the FIRES surveys iFranx et'aP 
|2P03; Forster-Schreiber et al. 2004), although there is a larger 
fraction of DRGs with z < 2 in the GOODS-S field compared 
to the FIRES samples (see §|5}. We believe the lower median 
redshift of the GOODS-S DRGs primarily results from the 
larger areal coverage combined with the somewhat brighter 
flux limit of the GOODS-S field. In addition, cosmic vari- 
ance between the GOODS-S and FIRES fields may contribute 
to this difference. 

However, roughly 20% (30/152) of the DRGs have redshifts 
z < 1 .5 (including two with spectroscopic redshifts), signifi- 
cantly lower than the typical redshifts reported for the FIRES 
samples.'^ We have inspected the SEDs of this low-redshift 

*^ lvan Dokkum et^I] 120031) find that one of the six DRGs in their spec- 



DRG subsample and found that the ACS through IRAC col- 
ors supports the derived redshifts with high fidelity for the 
majority (27/30) of galaxies. For these DRGs, the IRAC pho- 
tometry shows a fairly robust turnover at rest-frame 1.6 /im 
near the expected peak in the stellar emission, and the SEDs 
at the inferred redshift are otherwise consistent with the ACS 
and ISAAC colors. In the remaining three cases, their ACS 
photometry has low S/N, leading to dubious results. Ignoring 
these three galaxies has no effect on the analysis and conclu- 
sions derived here, but we include them in the sample for com- 
pleteness. In §|6l we conclude that the lower redshift (z < 2) 
DRGs are mostly heavily extincted starbursts, and are proba- 
bly part of the class of dusty EROs at z > 1, which typically 
have red J-Ks colors that satisfy the DRG selection criterion 
JSmail et alJ2002HFranx et al.l2003l) . 

2.2. Spitzer 24 jjLm Source Detection and Cataloging 

f apovich et al. (2004b) describe the data reduction 
and point-source photometry methods applied to the 
SpitzerMWS 24 image. They show that the GTO 
24 /im data reach a 50% completeness limit at 60 /iJy. We 
reanalyzed the simulations discussed in Papovich et al. to 
estimate the flux uncertainties as a function of MIPS 24 /im 
flux density. The median 24 jim S/N is Ri4 for objects with 
/i,(24/im) = 80-90 /iJy (the 80% completeness Hmit). We 
cross-correlated 24 /im sources with /i,(24/im) > 50 /iJy 
(the S/N w 3 limit) to the /T^-band catalog and identified 
matches with a radius of 2". Roughly one-half (74/153) 
of the DRGS are detected by SpitzerMWS at 24 /im. The 
majority of these 24 /im-detected DRGs (71/74) have IRAC 
counterparts associated with each DRG in the K^-hwd 
image. The remaining three appear associated with IRAC 
sources at distances 0."5 <r< 2"0 from the K^-hwd source. 
For the source density at the flux limit of the 24 /im data 
(Papovich et al. 2004b), we expect a random-association 
probability of 0.03 within a 2" radius, which is consistent 
with the three sources with no IRAC counterpart. Therefore, 
these sources are likely change alignments. 

Figure n shows the redshift distribution of the DRGs de- 
tected at 24 /im. The distribution of this sub-population is 
similar to those DRGs not detected at 24 /im, except for an 
apparent spike in the redshift distribution at z — 2.2 for the 
24 /im-undetected DRGs. The 24 /im-detected DRGs have 
a mean redshift (z) = 2.0, slightly lower than the mean red- 
shift for the DRGs with no 24 /im detection, (z) =2.3. Using 
a Kolmogoroff-Smirnov statistic, there is a moderate likeli- 
hood (90% confidence) that the two sub-populations of DRGs 
are drawn from different parent samples. Different redshift 
distributions in part may arise because the 24 /im subsam- 
ple includes a large number of lower-redshift galaxies. How- 
ever, at z ^ 2 the 7.7 /im emission feature from polycyclic 
aromatic hydrocarbons (PAHs) lies in the 24 /im bandpass, 
and thus one might expect more IR-detected galaxies to ap- 
pear. The sharp increase in the number of 24 /im-undetected 
DRGs al z — 2.2 may result from the fact that at this red- 
shift the DRG selection begins to pick up galaxies with strong 
Balmer/4000 A breaks that shift between the J and /C^j-bands. 
Such galaxies presumably have low specific SFRs (SFR per 
unit stellar mass), and perhaps have lower IR emission, so 
fewer would be detected at 24 /im. 

troscopic sample has Zspcc = 1.19. They estimate that the red J-Ks color 
selection has a ~ 20% contamination of (dusty) galaxies at lower redshift, 
probably consistent with the fraction of DRGs atz < 1.5 in our sample. 
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2.3. X-ray Source Detection and Cataloging 

The CDF-S has aim-point flux limits (S/N= 3) in the 0.5- 
2.0 keV and 2-8 keV bands of « 2.5 x 10"'^ erg cm'^ s"' 
and « 1.4 X 10"'^ erg cm"^ s"', respectively. The com- 
pleteness Umit over « 90% of GOODS-S field is « 
1.3 X lO-i*" erg cm'^ s'^ (0.5-2.0 keV) and « 8.9 x 
10"'^ erg cm"^ s~' (2-8 keV). Assuming an X-ray spectral 
slope of r = 2.0, a source detected with a flux of « 10"'^ erg 
cm"^ s"' would have both observed and rest frame luminosi- 
ties of « 5.2 X lO"*' erg s"' and w 7.6 x lO'*^ gj-g s"' at z = 1 
and z = 3, respectively, assuming no Galactic absorption. 

In this st udy we use the main an d supplementary Chandra 
catalogs of Ale xander et al. (2003). The median positional 
accuracy for the sources in the GOODS-S field in the main 
Chandra catalog is Q'.'6. We matched all ^Tj-band sources 
to the X-ray catalog within a radius of 1" and identified all 
matches in both the soft- and hard-band catalogs. Nearly 
one-seventh (22/153) of the DRGs have X-ray detections in 
either the soft or hard bands (or both). Of these, 12 are de- 
tected at 24 /im, while the remaining 10 are not. All of the 
X-ray detected DRGs have IRAC counterparts. 

2.4. High-Redshift Sample of IR-luminous Galaxies from the 

HDF-N 

We construct a comparison sample of galaxies within 
flie northern Hubble Deep Field (HDF-N, 'Willia ms etaP 
which spans all types of galaxies (not just the red- 
dest galaxies identified by the DRG-selection), and ex- 
tends our analysis to fainter 24 fim fluxes. This allows 
us to study how the DRGs are drawn from the general 
galaxy population at similar redshifts. Galaxies were se- 
lected from the NICMOS HDF catalog of M. Dickinson et 
al. (Dickinson et al. 2000; Papovich, Dickinson, & Ferguson 

I2OOL .Dickinson et al 2003) . and matched to the deep IRAC 

and MIPS 24 fj,m observations of the northern GOODS 
(GOODS-N) field (M. Dickinson et al., in preparation; 
R. Chary et al., in preparation). From the NICMOS-selected 
catalog, we identify 24 galaxies in the HDF-N with Ks < 
23.2 mag and redshift, 1.5 < z < 3.5, roughly in the same 
range as in the GOODS-S DRGs. Spectroscopic redshifts 
are available for 13 of the galaxies in this sample (see 
[Dickinson et al. 2003, and references therein). For the re- 
maini ng 11 galaxies, we u se the photometric redshift catalog 
from Budavari et alJ (l2000i) . Of this sample, 19 are detected 
with the MIPS 24 pm imaging to f^(24fj.m) > 10 /iJy. 

Figure Q] shows the redshift distribution of the HDF-N 
galaxies. We note that four of these galaxies with K, < 
23.2 mag have (/iio-^s)ab > 1-6, and possibly satisfy the 
DRG color criterion (allowing for differences in the HDF-N 
and GOODS-S filter sets). Of these, three have f^(2Apm) > 
50 /iJy, and so could have been detected in the 24 /im of the 
CDF-S. 

2.5. Low-Redshift Galaxy Samples from COMBO-17 

In §|6l we compare the SFRs and stellar masses of the high- 
redshift DRGs to those of lower-redshift samples. The GTO 
SpitzerlMWS 24 /im imaging intersects roughly 700 arcmin^ 
of the COMBO-17 survey (Wolf et al. 2003) in a substantially 
larger region encompassing the 130 arcmin- GOODS-S field. 
COMBO-17 provides phot ometric re dshifts for galaxies with 
R < 23.5 mag to z < 1.3 ( Wolf et al. 2004). Where possible, 
we replaced many of these photometric redshifts with sp ectro- 
scopic ones (Le Fevre et alJl2004l l\^zella et alJl2005h . We 



then constructed samples of galaxies with < 23.5 mag in 
two redshift sHces, 0.3 < z < 0.5 and 0.65 < z < 0.75. We 
cross-correlate the COMBO-17 catalogs with the GTO MIPS 
24 pm catalogs to identify matches within a 2" radius. We 
find 1495 galaxies in the full z ^ 0.7 sample, of which 464 
are detected in the MIPS 24 /im catalog to 24 /tm 50% com- 
pleteness flux hmit. Similarly, there are 1269 sources in 
the full z ^ 0.4 sample, with 276 sources detected in the 
MIPS 24 pm catalog (see also the discussion in Bell et al. 
2005; Le Floc'h et al. 2005). Each sample contains approxi- 
mately the same co-moving volume in these redshift intervals 
(~ 10^ hjQ Mpc^), although we find relatively fewer objects at 
z ~ 0.4 than z ^ 0.7 to the same magnitude limit. This results 
from large-scale cluste ring in this field, which is known to be 
underdense at z 0.4 (Wolf et al. 2003). While IRAC imag- 
ing also exists for most of the COMBO-17 field, we have not 
included it in our analysis of the galaxy stellar masses as it 
provides observations longward of rest-frame 2 pm , past the 
peak of the stellar emission (see discussion in § I6.4> . 

3. THE REST-FRAME OPTICAL AND NEAR-IR COLORS OF DRGS 

Galaxies at z ~ 2-3.5 with a strong 4000 A/Balmer break 
should hav e colors that s atisfy the (/-^Ovega > 2.3 mag cri- 
terion (Fra nx et "311120031) . This color selection is also sensi- 
tive to starburst galaxies at z > 1 whose light is heavily ob- 
scured by dust. Forster-Schreiber et al. ( 2004) note evidence 
for both galaxy types in the /gu-y,, Js-H, and H-Ks color 
distributions of DRGs, which is similar to the reputed nature 
of Bz/T-selected objects Daddi et al. (2004). We find that the 
full rest-frame UV to near-IR colors of DRGs provide further 
support for this dual population. 

In Figure|2lwe compare the SEDs of the GOODS-S DRGs 
with empirical and theoretical UV-IR SEDs for local galax- 
ies. The DRGs have colors that fit within the envelope de- 
fined by various galaxy types. Many of the DRGs have col- 
ors close to the locus of that for early-type galaxies. There 
is also a class of DRG with optical co lors redder than even 
the elHptical template of lColeman et aljfl980 ). Most of them 
are 24 /im sources, and they have very red UV-optical rest- 
frame colors, consistent with dust-obscured starbursts. The 
flux ratio between the 24 /im and near-IR emission is consis- 
tent with that of galaxies with massive starbursts, such as the 
local ULIRG Arp 220. Most of the DRGs show a prominent 
inflection in their SEDs at rest-frame 1 .6 pm as expected in 
the SEDs of galaxies dominated by the light of composite stel- 
lar po pulations (e.g., Simpson & Eisenhardt 1999; Sawicl^ 
l2002h . The ensemble DRG photometry in the mid-IR shows 
slight evidence for a peak at ^ 8 /im rest-frame coincident 
with the PAH 7.7 pm emission feature, suggesting that the 
emission from the majority of DRGs stems from stellar pro- 
cesses rather than nuclear activity. There is a slight decline 
in the flux density at 9- 10 /im, possibly consistent with sil- 
icate absorption, which is observed in both starburst galaxies 
and AGN. A small subset of the DRGs have red colors around 
A ~ 1 - 5 /im rest-fr ame, a sign ature of AGN emission (e.g ., 
lRiekell 978: Neugeba uer et all 1979. and discussion in § l6.3t . 

In Figure |3] we show an optical and near-IR color-color 
diagram to study the stellar populations and dust extinction 
in the DRGs. The simple models plotted in the figure bound 
the range of colors observed in the GOODS-S DRGs, com- 
parable to the findings of Labbe et al. ( 2005| for a smaller 
FIRES sample. DRGs with the blue Kg - [4.5pm] colors re- 
quire a substantial population of mature stellar populations 
with a strong 4000 A/Balmer break, which produces the 
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Fig. 2. — Ensemble photometry of the DRGs in the GOODS-S sample. The data points are the measured flux densities from the ACS (£435, Vgo6, '775, Z85o)> 
ISAAC (J, H, Ks), IRAC ([3.6/im], [4.5/xm], [5.8/jm], [8.0/xm]), and MIPS 24 /im images, normahzed to a common flux density at 1.6 fim rest-frame. Open 
symbols show the photometry for galaxies detected at 24 /im with MIPS, and solid symbols show galaxies with no 24 /im detection. Stars correspond to galaxies 
detected in the X-ray data. The shaded region is defined by the SEDs for ordinary galaxies from Coleman, Weedman, & Wu (1980; CWW80), and span Hubble 
types from elliptical galaxies through Sab, Sbc, and Scd spirals, to the bluest Magellanic irregular, Im, templates. The solid cyan curves show the SEDs from 0.1 
to 3 /im for stellar populations from the Bruzual & Chariot (2003; BC03) models formed with constant star formation for 100 Myr with dust extinction of 1.5, 
4.5, and 7.5 mag at 1600 A using the jCalzetti et al. i2000) extinction law. The solid b lue curve shows th e SED of the local ULIRG Arp 220 from an empirical 
model for 1-5 fim <Silva et alll998ft and a composite spectrum covering 5 — 3000 fim 1 Spo on et alJ2004i) . The latter shows the prominent PAH emission feature 
at 7.7 /im, and silicate absorption at 9.7 /im. The short-das hed, red curve show s the SED of the local ULIRG and AGN Mrk 231, and the long-dashed magenta 
curve shows that for the local Seyfert 2 galaxy NGC 1068 ILe Floc'h et all2(X)li) . both of which rise smoothly through the near- and mid-IR. 



red J—Ks color. DRGs with redder K^ — [4.5/im] colors re- 
quire ongoing star formation with substantial dust extinction. 
Changes in the redshift of the model stellar populations from 
z = 2.2 to 3.5 or 1.5 shift the expected colors by ±0.3 mag 
in Ks - [4.5^m] mag, respectively, but have little effect on the 
'775 -Ks colors. The models bound the majority of the colors 
of the DRGs (with the exception of several X-ray sources, 
whose optical to near-IR colors are possibly influenced by 
AGN). Thus, the colors of the DRGs imply a mix of old and 
young stellar populations, in some cases with substantial ex- 
tinction. The MIPS lAiiva data support this interpretation. 
The top panel of Figure|3]shows the ^'^-[4.5/im] colors of the 
DRGs detected at 24 /im and those with no 24 /im detection. 
The DRGs with 24 /im detections have redder Kg - [4.5/im] 
colors, suggesting that the emission from many DRGs is at- 
tenuated by a large dust opacity. 

Of the sample of 132 DRGs with IRAC detections, there 
are only 6-9 DRGs (within the photometric uncertainties) 
that have ijjs-Ks > 4 mag and /Tj- [4.5/im] < 1.5 mag (in- 
cluding two sources with 24 /im detections and one X-ray 
source). These color thresholds should identify passively 



evolving galaxies with ages greater than w 0.75 Gyr over the 
redshift range z ~ 2-3, and these colors bound the "red and 
dead" DRGs in the FIRES sample from the southern Hub- 
ble Deep Field (HDF-S) ( Labbe et al. 2005). Assuming the 
GOODS-S DRGs with these colors are uniformly distributed 
over the GOODS-S area and redshift range 2 < z < 3 implies 
they have a number density of 0.7 -2 x lO'^/i^g Mpc"^. This 
contrasts with ILabbe et alJ J2005ft who found a higher num- 
ber density, 1.9 x IQ-'^hjl Mpc-^ for DRGs in the HDF-S 
assuming the same redshift interval. The difference between 
the DRGs in these fields is not removed by relaxing the color- 
selection criteria further Based on the full photometry of 
our DRG sample (see §|5ll, we find that 15 DRGs (- 10%) 
have best fit models consistent with old ages (> 1 Gyr), lit- 
tle dust (E[B-V] < 0.1), and passive evolution. They have 
/775 -Ks > 2.6 mag and Ks- [4.5/fm] < 1.1 mag, both some- 
what bluer than the original limits. These galaxies suggest 
that the number density of massive, passively evolving galax- 
ies at 2 < z < 3 is 3.2 x lO'^/i^Q Mpc""*, still nearly an order of 
magnitude lower than that in the HDF-S. However, some of 
the discrepancy likely arises from the fact that the FIRES data 
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Fig. 3. — The Ks — [4.5fim] versus (775 —Ks, color-color diagram for the GOODS-S DRGs. Symbols show the derived colors from the ACS, ISAAC, and IRAC 
photometry. Black, open circles show DRGs not detected at 24 /im with MIPS. Red, open squares denote DRGs with 24 fim detections. Filled stars con'espond 
to sources detected in the soft or hard X-ray bands with Chandra. Arrows denote upper Umits on the ACS (775 — Ks colors. The lines show the expected colors 
of stellar populations at z = 2.2 (the median redshift in the DRG sample) using the Bruzual & Chariot |2003) models. The dashed line shows the colors of a 
passively evolving stellar population formed in a single burst with an age of 50 Myr to 2.8 Gyr; hash marks indicate the colors at ages of 0.5, 1, and 2 Gyr, as 
labeled. The solid lines show the colors of a stellar population forming with constant star formation for the same ages as above, but with a dust extinction of 2 
and 6 mag at 1600 A, as labeled. The arrow shows the expected differential color for 2 magnitudes of extinction at 1600 A. The histograms on the top and left 
panels show the color distributions of those DRGs detected at 24 /^m (red solid lines) and undetected at 24 ^m (shaded regions). 



achieve fainter near-IR flux densities, and parenthetically we 
note that two of the three candidates for passively evolving 
DRGs in the Labbe et al. sample have > 23.2 mag. Nev- 
ertheless, the difference supports the notion that these sources 
are highly clustered (Da ddi et al. 2003), and that the HDF-S 
itself has an unusual overdensity of them. The low number 
density of massive, passively evolving galaxies in our sam- 
ple supports the assertion that the densit y of passive galax- 
ies is rising strongly at z < 2 (see, e.g., iDaddi et al.ll2005at 
iLab be et al . 2005). 

There are few DRGs with /wg - < 2 (see also, 
ILabbe et al.ll20()5l) . Four of the GOODS-S DRGs have /vvs - 
Ks < 2 and Ks - [4.5 iim] < 0. One of these galaxies is an 
X-ray source at z ~ 2.3 whose colors may be affected by 



an AGN. The remaining three have redshifts z > 3 and rela- 
tively flat SEDs from UV to near-IR rest-frame wavelengths. 
The red J-Ks color arises from a weaker 4000 A/Balmer 
break apparently augmented by photometric errors or emis- 
sion lines in the passbands,'^ which push the /-band fainter 
and Ks-hand brighter by small amounts, and conspire to pro- 
duce (J-Ks)Yeg!i > 2.3 mag. In general we retain these sources 
in the sample for completeness (although they have a negligi- 
ble effect on our results). In our analysis below, we also con- 
sider a restricted sample with 1.5 < z < 3.0, thus excluding 

" van Dokkum et al. (2004) find small emission-line corrections of ~ 
0.1—0.2 mag to the ^fj-magnitudes in a sample of seven, bright (^^[Vega] < 
20 mag) DRGs. 
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these objects. 

4. TOTAL LUMINOSITIES AND STAR-FORMATION RATES OF 
HIGH-REDSHIFT GALAXIES 

4.1. Estimating the Total Infrared Luminosities of 
Star-forming Galaxies 

At the redshifts of the DRGs, z 1 -3.5, the Spitzer 24 fim 
probes the rest-frame mid-IR, which broadly correlates with 
the total thermal IR lumino sity, Lir = L(8 — lOOOA^m) (e.g., 
SpinogHo et al." "1995; Roussel et al. 2001"; "Charv & Elbaz 
2001: Dale et al.. 2Q01j Elbaz et al. 2002: Papoyich & Bd] 
2002) ■ We convert the observed 24 /j,m flux density to a rest- 
frame luminosity density at 24/(1 +z) ^m. We then correct 
these values to a total IR luminosity using the Dale & Helo3 
(|2002) IR template SEDs assuming that a given rest-frame IR 
luminosity density translates uniquely to a single SED tem- 
plate. If we instead used the IR templates of Charv & Elba3 
J2001h . then we would derive IR luminosities a factor of 
2-3 higher relative to those of Dale & Helou for galaxies 
at LiR ~ lo'^^"'^ Lq (with smaller differences for lower- 
luminosity galaxies). A recent study of IR-luminous galaxies 
at z < 1.2 indicates that IR-luminosities estimated from the 
Chary & Elbaz models have a scatter of a factor of two com- 
pared t o IR luminosities der ived from the radio-far-IR corre- 
lation (Marcillac et al."2005'). In contrast, the IR luminosities 
estimated using the Dale & Helou models provide a tighter 
correlation with IR-luminosities derived from the radio-far- 
IR correlation, with a scatter of 40%, suggesting these tem- 
plates possibly better reflect reality. Some scatter is inherent 
in this estimation of the total IR luminosity: Chapm an et al.l 
(|2003) find that the temperature-luminosity distribution in 
IR-luminous galaxies has a scatter of roughly a factor of 2- 
3 in IR lu minosity fo r galaxies with fixed dust temperature. 
Howe ver. IPaddi et al . (2005b) find that the (Chary &Elba3 
12001'^ IR model template with Lir = 10'^-^ L© fits the average 
SED of 24 /im-detected BzK objects at (z) = 1.9, suggesting 
that the uncertainty in the templates is not severe. Neverthe- 
less, we add 0.3 dex as a systematic error on the inferred Lir 
to account for the systematic scatter in this conversion. 

The uncertainty of the photometric redshifts leads to 
another source of uncertainty in the conversion from 
Liy(24/im/[l to the total IR luminosity. Owing to the 
large bolometric corrections from the mid-IR to the total IR 
luminositie s, small changes in the redshift have a significant 
effect (see Pap ovich & Belli 120021) . We find that taking the 
68% confidence range on the photometric redshifts of the 
DRGs leads to variations in the inferred Lir of 0.4 dex. We 
add this source of error in quadrature with the uncertainty 
from the IR templates, bringing the total error budget on the 
derived IR luminosities to 0.5 dex. 

Figurel^shows the total IR luminosities of the DRGs. The 
completeness limit (converted to a total IR luminosity in the 
same way as for the 24 /im detected galaxies) is indicated in 
Figure 0] as the solid line. The figure also shows the total IR 
luminosity for the HDF-N galaxies in this re dshif t range de- 
rived from deeper 24 /im observations (see § I2.4> . The total 
IR luminosities of the 24 /im-detected DRGs are lO""'"* Lq, 
or 10""'^ ^' Lq if we exclude X-ray detected sources. The 
majority of these objects have IR luminosities comparable to 
local ULIRGs, Lir > 10'^ Lq, which if attributed to star- 
formation implies SFRs greater than 100 A4q yr"'. 

Several of the DRGs (6/152) have IR luminosities, Lir > 
lO'^ Lq (so-called Hyper luminous IR galaxies, HyLIRGs). 
These IR luminosities are comparable to those of PG quasars 
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Fig. 4. — Total IR luminosities, Lir = L(8 — lOOO/xm), of galaxies inferred 
from their observed MIPS 24 ^m emission. Red circles show the IR lumi- 
nosities for the GOOD-S DRGs detected at 24 /tm, black stars show DRGs 
with both X-ray and 24 fim detections, and blue squares show objects de- 
tected by MIPS 24 fim in the HDF-N. The solid line denotes the 50% com- 
pleteness limit of the GTO 24 fim data in the CDF-S, f\y(2Afim) = 60 /jjy. 
Error bars include only uncertainties on the 24 fim flux density. We estimate 
that the systematic uncertainties on the conversion between /y(24/im) and 
L(8- 1000/im) are typically RJ 0.5 dex (inset, and see text). 



at z > 1, which have warm thermal dust temperatures 
JHaas et al.ll2003l) . Locally, HyLIRGs — and many ULIRGs 
with Lir > lO'^ -^ Lq have rest-frame optical emission spec- 
tra characteristic of Seyfert galaxies (e.g. Soiferet al. 19951 
I Veilleux etai] fr995: Sanders & Mirabel 1996; Veifleux et aQ 
1 1999!), and the IR emission possibly originates from AGN 
processes. We suspect that the most IR-luminous DRGs 
may have a contribution to their bolometric emission from 
AGN. The majority (5/6) of the DRGs with infeiTed Lir > 
10'"^ Lq are detected by Chandra, compared to the 10% X- 
ray-detection fraction over the whole sample. This X-ray de- 
tection fraction is consistent with the limit on the AGN frac - 
tion for the coeval sub-mm galaxies (I Alexander et alJl2005l) . 
a source population that also has inferred Lir ~ lO'^ Lq. 
HyLIRGs are also present in high-redshift IRAC-selected 
AGN samples, a nd ^ 50% ar e undetected in deep X-ray data 
( Alonso-Herrer o et alJl200.5l) . Therefore, the X-ray emission 
from AGN in many of these IR-luminous objects may be 
attenuated by dust below the detection limit of the surveys. 
However, if AGN contribute to the emission in DRGs with 
^'iR ^ lO''' Lq, then we may be overestimating the galaxies' 
IR luminosity. Although the Dale & Helou (2002) IR tem- 
plates include galaxies with Lir > lO'-' Lq, using a template 
for Mrk 23 1 with a known AGN and warmer dust temperature 
would reduce the inferred IR luminosity by factors of ~2-3. 
To limit the effects of any bias caused by IR template uncer- 
tainties for the highest luminosity DRGs, we consider below 
how restricting our sample to galaxies without X-ray detec- 
tions, Lir < lO'^* Lq, and IR colors indicative of AGN (see 
§ I6.3> affects our analysis. 

4.2. The Relation Between the UV Spectral Slope and the 
Infrared Excess 

Locally, UV-selected starburst galaxies have a relation be- 
tween their UV spectral slope, /3, where f\ ^ X^, and the 
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ratio of their IR to UV luminosit y (also termed the "infrared 
excess '; e.g.. lMeurer eTallfTQQQl) . This relation links an ob- 
served increase in (3 with greater dust extinction of the intrin- 
sic spectrum of a young stellar population. The dust-absorbed 
UV light is radiated in the thermal infrared, yielding an anti- 
correlation between the UV and IR. While this relation holds 
for a range of UV-luminous starburst galaxies jMeurer et alJ 
Il999h . it does not apply to other types of star-forming galax- 
ies, in cluding normal galaxies (Kona et al. 2004; Buat et al. 
I^POS), and ULIRGs (Goldader et al. 2002), apparently due to 
geometry effects between the U V an d IR emitting regions. 
Using the IR luminosities from § 14.11 we test whether the lo- 
cal relation between the UV spectral slope and the IR excess 
applies to the DRGs. We also test if it applies to other high- 
redshift galaxies detected by MIPS at 24 fim. 

Following Meureretal. (1999), we calibrated measure- 
ments of l3 by comparing the UV spectral slope derived from 
the full UV SEDs of local starburst galaxies, /3spec, with that 
measured using only a broad-band photometric color, /3phot- 
The correction for the difference between P^psc and /?phot is 
necessary as starburst-galaxy spectra contain numerous ab- 
sorption and emission features. A spectroscopi c UV slope is 
fit to spectral windows to avoid these features iCalzetti et al.l 
[1994), but they are unavoidable in measuring (3 from broad- 
band colors. We fit the spectral slope, /3sDec, using the 
UV SEDs of starburst galaxy templates (Kinney et al. 19961, 
which have dust extinction varying from E{B-V) < 0.10 to 
0.61 < E(B-V) < 0.70. We then shifted these spectra to the 
observed frame at z = 0.5 -4, and measured the observed ACS 
^435 -^606, B435 -'775, ^606-/775, and WFPC2 y606-/8i4 Col- 
ors. We measured the approximate UV spectral slope from the 
ACS and WFPC2 broadband photometry, /Sphoi, using the ef- 
fective wavelengths of the filters, /3phot = 3.20(V606-^8i4)-2.0 
for the WFPC2 bands, and /3phot = 2.91(B435 - Veoe) - 2.0, 
/3phot = 1.59(B435-;775)-2.0, and /3phot = 2.14(y606-Z85o)-2.0 
for the ACS bands. We then fit a quadratic polynomial to the 
difference, /3spec - Pphot, as a function of redshift. 

To derive the UV spectral slopes of the DRGs, we want 
to span the longest wavelength baseline in the rest-frame 
wavelength range 1250-2800 A. At wavelengths longer than 
2800 A there may be a significant contribution to the SED 
from A- and later-type stars from previous star-formation 
episodes (e.g., I Calzetti et al.lll994l) . Therefore, for galaxies 
at z > 2.2 in the GOODS-S field, we will use the Veoe-^sso 
color to measure f3. For galaxies with 1.7 < z < 2.2, we will 
instead use the B435-/775, where the lower redshift bound 
results from the fact that below this redshift the /775-band 
probes A > 2800 A. For galaxies with 1.2 < z < 1.7, we 
use B435 - Veoe to derive the UV spectral slope. For z < 1.2, 
there are no ACS colors that measure the UV spectral slope at 
A < 2800 A. 

For galaxies from the HDF-N with WFPC2 photometry 
and redshift, z, we derive the UV spectral slope, /?, using the 
quadr atic fit to the empirical relationship (cf., .Meurer et al. 

P = 3 .20(Vm - /814) - 4.45 + 2 .03z - 0.423zl ( 1 ) 
For the DRGs with ACS photometry we use 

( 2.91(B435-V6()6)-3.35+1.49z-0.438z^ 1.2<z<1.7 
P=< 1. 59(^435 -/775)- 3. 92 +1.97Z- 0.5 12z^ 1.7<z<2.2 
[ 2.14(y606-Z85o)-4.76 + 2.14z-0.418z2, z > 2.2. 

(2) 



We find that the different formulae give consistent results 
within the photometric errors for galaxies in redshift intervals 
where several ACS colors can be used to measure Using 
these formulae, we then calculate the UV spectral slopes of 
the DRGs and HDF-N galaxies using the appropriate ACS 
and WFPC2 bands, respectively. We also use uncertainties 
of the observed colors and the covariant uncertainties of the 
polynomial fits to estimate the error in f3. 

Figure |5] shows the relation between (3 and the IR/UV lu- 
minosity ratio for the high-redshift galaxies, and compares 
them to the local relation from Meurer et aL (1999). Many of 
the HDF-N 24 /zm-detected galaxies at 1.5 < z < 3.5 have 
UV spectral slopes and IR/UV luminosity ratios that lie near 
the local relation. These systems have IR luminosities in the 
range Lir ^ 10'*^"" Lq, comparable to the IR luminosities in 
the Meurer et al. sample. There is a trend for galaxies with 
higher IR/UV luminosity ratios to move away from the lo- 
cal relation in the sense that they have more IR luminosity 
than otherwise predicted from their UV luminosity and spec- 
tral slopes. 

This trend in galaxies with higher IR luminosities hav- 
ing higher UV/IR ratios is continued for the 24 /im-detected 
DRGs. The DRGs span /3 ~ -2 to 2 (with the exception of one 
object with /3 ~ 3.7 and large photometric uncertainty). Many 
of DRGs with Lir < 10'^ Lir and some with Lir > 10 Lq 
have IR/UV ratios near the Meurer et al. relation. However, 
most of the DRGs have IR/UV luminosity ratios of more than 
one order of magnitude in excess of what would be predicted 
from their spectral slopes. That is, the total amount of star 
formation in these galaxies will be underestimated from their 
UV rest-frame luminosities and spectral slopes alone. This 
result is qualitatively unchanged if we restrict the DRG sam- 
ple to redshifts 1.5 < z < 2.5, where the total correction fac- 
tor from observed f^(24^m) to Lir spans a relatively narrow 
range compared to that for the full sample. 

The DRGs occupy a very similar range of (3 and IR / UV flux 
ratios as has been observed for local ULIRGs (Goldade r et alJ 
l2002h . In local ULIRGs the geometry of the forming star clus- 
ters and dust is highly complex. The UV and IR emitting re- 
gions are typically displaced from one another, or "patchy", 
such that the regions that dominate the UV emission are unas- 
sociated with the regions producing the large IR luminosity 
(as is the case for local IR-luminous galaxies with ISO IR and 
HST UV imaging; see ICharmandar is et alJl20()4l) . A similar 
situation probably holds for the 24 /im-detected DRG pop- 
ulation. It is also plausible that the galaxies with the largest 
IR luminosities have more complicated geometries to account 
for the fact that these galaxies have the largest offsets from 
the local relation. Charmandaris et al. noted that for many of 
the local LIRGs, the ISO source is offset from the regions that 
dominate the UV emission by as much as several kpc, and 
such objects can have large IR excesses. This effect should be 
even stronger for ULIRGs. At redshifts typical of the DRGs, 
this corresponds to less than 1", and will be unresolvable with 
MIPS at 24 /im. 

Most of the X-ray-detected DRGs have IR/UV luminosity 
ratios and UV spectral slopes comparable to the rest of the 
DRG sample, although several have the most extreme IR/UV 
luminosity ratios or UV spectral slopes (logLiR/L[1600A] ^ 
4, and/or (3 ~ -2), lying away from the other DRGs. We sus- 
pect that an AGN contributes substantially to the I R emission, 
the UV emission, or both. In this case, the iDale & Heloul 
( I2002t IR templates may overestimate the IR luminosity in 
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Fig. 5. — Relation between the IR/UV luminosity ratio and the UV spectral slope, /3, where f\ ^ A''. The solid line shows the relationship for local 
UV-luminous starburst galaxies from Meurer et al. 1 1999). Red pentagons denote the GOODS-S DRGs, black stars indicate DRGs with X-ray detections, and 
the symbol size scales with IR luminosity (see plot inset). Blue squares denote galaxies with 24 fj,m detections from the HDF-N. The top abscissa shows the 
extinction at 1600 A that corresponds to /3 (see Meurer et al. 1999). Error bars on the IRAJV luminosity ratios only correspond to the uncertainties on flux density 
measurements. The systematic uncertainty in the conversion between /^(24/im) and Lir is RJ 0.5 dex, and is indicated by the inset error bar. 



these objects (see § 14. 1> . 

4.3. Inferring the Star Formation Rates of High-Redshift 
Galaxies from Bolometric Luminosities 

Nearly all of the bolometric luminosity from star-forming 
regions is emitted in the UV and IR (e.g., Bell 2003j). 
Therefore, we estimate the instantaneous star-formation rates 
(SFRs) for galaxies in our samples using the combination of 
their UV and IR luminosities. 

We use the SFR conversion from lBell et al.l (l2QQ^dl. base d 
on the UV and IR calibration presented bv lKennici 

^/Mq yr-i = 1.8 X 10"'° x (Lir + 3.3 Lzsoo)/^©, (3) 

where Lir is the total IR luminosity, and L280() = v L^(2%QQA.) 
is the monochromatic luminosity at rest-frame 2800 A inter- 
polated from the ACS photometry. We have adjusted the SFR 



to correspond to a single power-law, Salpeter IMF with mass 
cutoffs of 0.1 and 100 7\/(0. Bell (2003) tested these UV+IR- 
derived SFRs against extinction-corrected Ha and radio- 
derived measures, finding excellent agreement with < 0.3 dex 
scatter and no offset. Nevertheless, it should be noted that 
the UV and IR-calibrations are based on local galaxy corre- 
lations. Although the indications are these hold at higher red- 
shifts (see, e.g., Elbaz et al. 2002 HAppleton et all2004i) . they 
should be used with some caution. 

This calibration explicitly assumes that star-formation pro- 
cesses account for the bolometric UV and IR emission. How- 
ever, the presence of AGN may also contribute some (or all) 
of this bolometric emission from accretion processes onto su- 
permassive black holes (SMBHs). Therefore, the SFRs de- 
rived using equation|3lmay be upper limits if AGN are present 
(modulo uncertainties in the conversion from the 24 /im to IR 
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Fig. 6. — Comparison between the SFRs derived from the extinction- 
corrected UV luminosity and the SFRs derived from the sum of the UV and 
IR luminosities as described in the text. Red circles denote DRGs detected 
at 24 /im, and gray triangles show upper limits for DRGs with undetected at 
24 fim. Black stars indicate DRGs detected in X-rays. Blue squares denote 
galaxies from the HDF-N at 1.5 < z < 3.5. Error bars on the UV + IR- 
derived SFR are ^ 0.5 dex, as indicated by the inset error bar. The diagonal 
lines indicate constant ratios of 1, 10, and 100, as labeled. 



luminosity). 

Figure |6l compares the SFRs derived from the sum of 
the UV and IR emission using equation 13 versus those de- 
rived using solely the UV rest-frame luminosity corrected 
for extinction usi ng the derived spectr al slopes (see e.g., 
iMeurer et al.lll999tlAdelberger & Steide]k.200 0). To calibrate 
the UV-derived SFRs, we compared the 1600 A luminos- 
ity (averaged ov er the range 1400-1800 A) to the SFR from 
iBruzual & Chariot 12003 ) spectral templates over a wide va- 
riety of star-formation histories. For a single power-law 
Salpeter IMF with mass limits 0.1 and 100 A4q, we find that 
a SFR of 1 A^Q yr"' corresponds to a luminosity density of 
Z,y(1600A) = 8.7 X 10^^ erg s"' Hz"^ for galaxies with ongo - 
ing star formation for > 10 Myr (see, e.g.. lMadau et alJl998h . 
Figure |6l shows this relation for the GOODS-S DRGs, and 
the HDF-N galaxies. The HDF-N galaxies extend to smaller 
SFRs, which are in closer agreement with the rest-frame UV 
data. The DRGs and many of the HDF-N galaxies with larger 
IR luminosities have significantly larger SFRs derived from 
the UV-i-IR than those estimated from the UV only. The dif- 
ference can be up to two orders of magnitude. 

5. STELLAR POPULATIONS AND STAR FORMATION IN HIGH 
REDSHIFT GALAXIES 

In § |3j we showed that the UV, optical, and near-IR col- 
ors of /-/Tj-selected galaxies are consistent with a multi- 
variate population of heavily dust-enshrouded starbursts and 
galaxies whose rest-frame optical and near-IR light are dom- 
inated by later-type stars. Here, we extend this analysis 
by comparing the full photometry of the DRG and HDF- 
N samples to stellar population synthesis models. As in 
[Papovich, Dickinson, & Ferguson (2001), we first consider a 
model of a single, monotonically evolving stellar population 
with a SFR that decays exponentially with a characteristic 
e-folding timescale, t. In reality, the star-formation his- 



tories of high-redshift galaxies are presumably more com- 
plex, involving stochastic events from mergers, interactions, 
feedback from star formation and AGN, as well as quies- 
cent star formation ( e.g., Somerville, Primack, & Fabet 200 J 
iNagamine et al. 2005; De Lucia etal. 2005). The models 
here should be considered as fiducial averages of past star- 
formation histories. Our definition of galaxy "age" is the 
time since the onset of star formation. Our monotoni- 
cally evolving models continuously produce new stars with 
young ages (albeit at a lower rate than in the past). For 
example, under our definition a stellar population formed 
with a constant SFR has age, f, while the mean age would 
be /5'(f)f(if/ /*(f)c/f = f/2, and the luminosity-weighted 
mean age (weighted heavily toward the short-lived early- 
type stars) would be younger still. The definition of galaxy 
age also neglects all previous discrete, episodes of star- 
formation. Older stellar populations from past star-forming 
events may very well exist, but be lost in the "glare" of 
the nascent stars. We consider these effects by using a 
second model that adds a maximally old stellar population 
formed in a single burst at z = cxo to the stellar popula- 
tions formed with the sim ple exponentially decaying models 
described above (see also. 'Pa povich. Dickinson. & Ferguson! 
2001 ; Dickinson et al. 2003). The latter model has a maximal 
stellar-mass-to-light ratio. We will use the single-component 
and two-component models to constrain the range of stellar 
masses and star-formation histories. 

5.1. Fitting the Models to the Photometry 

We fit the galaxy photometry with the Bruzual & CharloJ 
12QQ3.) stellar-population synthesis models. While 
the metalhcities of the DRGs are poorly known 
jForster-Schreiber et alJ 12004 . |van Dokkum et al.' ( 12004 
provide evidence for s olar and s uper-solar metallicities 
for luminous DRGs. 'Sha plev et al.i (i2004i 120051) have 
estimated the metallicities of massive (A4 > 10'^ Mq) 
[/„-dropout-selected LBGs at similar redshifts to be approx- 
imately solar. Because the expected colors for the majority 
of massive I/„-dropouts satisfy the D RG selection criteria 
("Redd v et alJ2005HShaplev et alJ20 05^. the solar-metallicity 
assumption for DRGs is also reasonable. Using different 
metallicities will affect our fitting results. However, the 
derived stellar masses vary by factors of less than 2-3 
(jPapovich. Dickinson. & Ferguson 2001). 

We use models with a single power-law Salpeter IMF with 
mass limits 0.1 and 100 Al©. Changing the shape of the IMF 
affects the derived stellar masses (and other stellar popula- 
tion parameters as well; see Papovich et al. 2001). For exam- 
ple, the Chabrier, Kennicutt, or Kroupa IMFs have a turnover 
in the mass function below 1 solar mass and produce stel- 
lar populations with roughly the same colors but with a stellar 
mass of 0.25 dex lower than that for the adopted Salpeter IMF. 
Although these other IMFs possibly better reflect nature, we 
choose to use a Salpeter IMF to facilitate comparisons of our 
results with those in the literature. As of yet there is no reason 
to expect the IMF to differ strongly from that observed in lo- 
cal galaxies (e.g. Baldrv & Glazebrook 2003; Larson 2005), 
although some empirical studies and theoretical predictions 
suggest a steepe r IMF may be required in high-redshif t mas- 
sive starbursts tFerguson et alJ2002t 'Baugh e t alJ2005h . 

We generate a^uite of photometry from the 
IBruzual & Chari ot (2003) models for galaxies spanning 
the range of redshifts in our sample with a redshift step- 
size of & = 0.05. The models range in stellar population 
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Fig. 7. — Illustration of .spectral synthesis model fitting results for one of the GOODS-S DRGs which has an indication of both young and old stellar 
populations. The star formation history used here is parameterized as a monotonia, decaying exponential with e-folding time r, as described in the text. The 
top panels show the 68 and 95% confidence intervals on various quantities plotted against the stellar population age: dust extinction (parameterized as the 
color excess, E[B — V]), and the stellar mass. The cross hairs show the most-likely parameter values in each two-dimensional projection of the full probability 
distribution function. The bottom panel shows the best-fit model spectrum in the observed frame, and the best-fit set of parameters for this galaxy. The data 
points show the ACS B435V606'775Z850, ISAAC JHKs, and IRAC 3.6-8.0 /xm photometry and errors. 



age from 10'' to 2 x lO'*' yr in the quasi-logarithmic 
steps prov ided in the m odels. We include dust extinction 
using the iCalzetti et al.l (1200 0) law with color-excess val- 
ues E{B-V) = 0.0-0.6 in steps of SE(B-V) = 0.025. 
For the Calzetti et al. extinction parameterization, these 
color excesses correspond to extinctions in the UV (rest- 
frame 1600 A) of A 1 600 = 0- 6 mag in increments of 
<5Ai6oo = 0.25 mag. We make the assumption that in these 
high-redshift galaxies the Calzetti et al. extinction law ap- 
plies to the stellar populations dominating the light observed 
in the ACS, ISAAC, and IRAC passbands (although it may 
not apply to th e ion izing source responsible for the far-lR 
emission, see § 14. 2> . Using different extinction laws would 



affect the derived stellar-population ages and extinction, 
but would not strong ly chan ge the infe rred stellar masses 
iPapovich. Dickinson. & Fergusoril 1200 Ih . We also make 
the assumption that the stellar populations dominating the 
rest-frame UV to near-IR also dominate the stellar mass. If 
a substantial fraction of the galaxies' stellar mass is obscured 
from view, then the stellar-mass derived from the SED 
modeling will be underestimated. We first allow for a range 
of star-formation histories with a SFR parameterized as a 
decaying exponential with a e-folding time, r, where the SFR 
at any age, t, is given by ^'(f) ^ expi-t/r). In our models, 
r ranges from 1 Myr to 100 Gyr in quasi-logarithmic steps. 
Short-duration t values correspond to instantaneous bursts 
of star formation while long-duration t values correspond 
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Fig. 8. — Same as Figure^ but for one of ttie DRGs witli a near-power-Iaw-Iike rest-frame UV-to-near-IR SED. Tlie best-fitting models generally prefer a 
SED that is dominated by young stars and heavily extincted by dust. 



approximately to constant star-formation histories. 

For each DRG, we fit all available bands among the 
ACS B435V606;775Z850, ISAAC JHK„ and IRAC [3.6^j.m] and 
[4.5/im] data. In the model fitting, we add a a/fi, = 4% error 
in quadrature to the photometric uncertainties on each band 
to account for the fact that the population-synthesis models 
do not continuously sample the rnodel p arameter space (see 
iPapovich. Dickinson. & Fergusonll2001h . The fits provide a 
normalization between the photometry and the model, and a 
minimum for each particular model with a distinct set of 
parameters. The IRAC [5.8/im] or [8.0/^m] photometry be- 
cause the longer-wavelength IRAC data generally have lower 
signal-to-noise ratios, and it is possible that at rest-frame 
> 2 fim PAH features or emission from obscured AGN shift 
into the IRAC bandpasses. Some galaxy SEDs show evidence 
for such features in these data (see Figure |2j, and so we ex- 
clude these points to avoid any potential bias. However, we 



find that a majority of the best-fit models broadly reproduce 
the [5.8/im] and [8.0/j,m] photometry, which lends credence 
to the fits. 

The model with the minimum chi-squared value, Xq, is the 
model with the best-fit for a given set of parameters (t, t, 
E[B-V], and A^; the latter is derived using the stellar-mass- 
to-light ratio of the model, see Papovich et al. 2001). Using 
the Ax^ difference between the chi-squared value derived for 
other models with a different set of parameters and Xq, we can 
construct confidence regions on each parameter in the model. 
For each galaxy we generated up to 1000 Monte Carlo real- 
izations for the photometry by perturbing the measured flux 
densities by a random value taken from a Normal distribu- 
tion with a standard deviation equal to the flux-density er- 
rors. We then refit the new photometry and re-obtain best fits 
on all model parameters. We do not take into account errors in 
the photometric redshifts, which can affect the derived star- 
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formation histories but generally have little effect on the in- 
ferred stellar masses (see Dickinson et al. 2003). We identify 
the Ax^ values from the fit to the measured data that encom- 
passes 68% and 95% of the best-fit values from the Monte 
Carlo realizations, which provides the equivalent confidence 
range on the model parameters for each source. 

5.2. Stellar population models with single exponentially 

decaying star-formation histories 

Figures [7] and |S] present examples of the results for two 
of the DRGs using the single-component star-formation his- 
tory models. The examples are typical of the two popula- 
tions identified by the (J-Ks)veoi > 2.3 mag selection. Each 
figure shows the probability distribution function in two- 
dimensional projections of the age-extinction, and age-mass 
planes. The figures also show the best-fit model spectrum 
over all parameters (i.e., the model with the minimum from 
all models) overplotted on the data. 

Figure0shows a galaxy with a strong 4000 A/Balmer break 
between the J- and /Tj-bands presumably due to a substantial 
mature stellar population. The galaxy also has a small amount 
of ongoing star-formation that produces the rest-frame UV 
light. The best-fit model corresponds to a star-formation his- 
tory that has undergone roughly f/r ~ 2 e-folding times. For 
this model, the galaxy formed most of its stellar mass in the 
past at a substantially higher SFR. The early-type (OB type) 
stars from the early onset of star formation have died off, and 
thus most of the stellar mass resides in later-type stars, which 
formed well in the past and now dominate the optical and 
near-IR rest-frame light. The best-fit model requires sub- 
stantial dust extinction: > 2 mag at 1600 A (68% confidence). 
The estimated stellar mass is quite robust for this set of star- 
formation histories. The 68% confidence range on the stellar 
mass ranges from 1.6-2.7 x 10" -Mq, with a most-likely 
value of 2.4 x lO" Mq. 

Figure|8lshows a galaxy with a heavily extincted recent star- 
burst. The best-fit model to this galaxy is a young stellar pop- 
ulation formed roughly instantaneously in a burst 40 Myr in 
the past. However, the 68% confidence range on the stellar- 
population age ranges from ^ 30 Myr to several Gyr. For 
all possible ages, the modeled star-formation history has un- 
dergone many e-folding times and has substantial dust ex- 
tinction, which is required to produce the red UV-to-near-lR 
rest-frame colors. This best-fit model requires 5 mag of ex- 
tinction at 1600 A. Owing primarily to the larger model de- 
generacies in the age and dust extinction for this galaxy, the 
stellar mass is less well constrained compared to the exam- 
ple in Figure The 68% confidence region on the stellar 
mass is 9 X lO** to 8.4 x 10'° Mq, with a most-likely value 
of 1.6 X 10'° Xq. 

5.3. Stellar population models with double component 

star-formation histories 

We also fit the DRG photometry models with a two- 
component star-formation history characterized by a pas- 
sively evolving stellar population formed in a previous "burst" 
with Zform = oo, summcd with the exponentially-decaying- 
SFR model above. For these models, the stellar-population 
age is the time since the onset of star formation in the mono- 
tonically evolving component. These models check the ef- 
fects of discrete bursts on the derived parameters. Our choice 
of a burst at Zform = cxd is a proxy for bursts at all times before 
the observed redshift, and reducing Zform would not strongly 



affect our conclusions. Placing the burst at Zform = oo provides 
a limit on the maximum stellar mass because it has the maxi- 
mal mass-to-light ratio possible at the observed redshift. 

Figure |9] shows the results for the two-component model 
fits to the two DRGs described above. For the galaxy shown 
in Figure a very young (f ^ 10 Myr), heavily extincted 
(^1600 ^ 3 mag) starburst dominates the UV rest-frame emis- 
sion. The optical and near-lR rest-frame emission is pro- 
duced by the previously formed population, which is now 
quite old. However, the 68% confidence range on the age 
and dust extinction shows that a large region of the param- 
eter space fits the data equally well. This is in contrast to 
the region permitted by the single-component model fits in 
Figure Clearly the interpretation of the ages and extinc- 
tion depend strongly on the assumed star-formation history. 
However, the derived stellar masses are fairly robust. The 
total stellar mass is nearly unchanged in this model relative 
to the single-component model above; the 68% confidence 
range is 1.6-2.3 x lO" M.q. The upper range of the confi- 
dence region drops simply because in the previous model we 
did not enforce the constraint that the age of the model be less 
than the age of the Universe at the given redshift. Such a con- 
straint is imposed (by construction) on the two-component 
fits, and as a result there is less time for stellar populations to 
evolve, produce large mass-to-light ratios, and increase the 
stellar mass. 

The best-fit two-component star-formation history model 
for the galaxy in Figure |8] shows little change compared to 
the single-component fit. In this case, young, reddened stars 
dominate the light at most wavelengths for both types of mod- 
els. Thus, for galaxies of this sort, adding the second compo- 
nent only provides an upper limit on the inferred stellar mass 
content. 

In general, the two-component fits for the DRGs appear 
more consistent with the data than the less complex, single- 
component models. Figure [lOlcompares the reduced Xq from 
the fits to the single-component models to the reduced Xo 
from the fits to the two-component models for all the DRGs. 
For many DRGs (91/152), adding the second model compo- 
nent to the fit has little effect on the minimum reduced x^, and 
these points lie near the unity relation indicated in the plot. 
For a large fraction of DRGs adding the second burst lowers 
the reduced minimum x^ significantly. It is possible that some 
of this effect results from the fact that we have not included 
the contribution of emission lines in the various bandpasses. 
However, as noted earlier, this contribution should be small, 
contributing < 0.1 ma g based on the observed emis sion- 
line equivalent widths (Ivan Dokkum et all 120041 12005b . ex- 
cept perhaps for galaxies with AGN, for which the models 
may not apply. Thus, in general this effect does not account 
for the smaller reduced x^ values observed in our galaxies. 
For example, the best-fitting single-component model for the 
galaxy in Figure deviates from the 7-band photometry at 
~ 1 jim by nearly 2 cr, whereas the two-component model is 
better able to reproduce this data point (see Figure|9j. 

The photometry in many of the DRGs is better represented 
by star-formation histories that are more complex than the 
simple, monotonically evolving exponentially decaying SFR. 
Formally, the x^ statistic rejects the single-component model 
for 13/152 DRGs at the 3cr level, while not rejecting the 
double-component model at this significance level. In con- 
trast, the x^ statistic never rejects the double-component 
model in favor of the single-component model at this sig- 
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Fig. 9. — Best-tit stellar-population models with two components to ttie two example DRGs shown in FigurelTlandlsl The solid grey line shows the best-fitting 
two-component model, and the data points correspond to the observed photometry. The two galaxy components correspond to a monotonically evolving stellar 
population (blue, short-dashed line) and a stellar population formed in a single burst at z = oo (red, long-dashed line). The inset best-fitting parameters refer to 
the monotonically evolving stellar population. The mass fraction is the percent contribution of the younger, monotonically evolving stellar population to the total 
stellar mass. The top panels show the 68 and 95% confidence regions on the age and dust extinction parameters and coiTespond to the galaxy in the spectrum that 
lies below each plot. The cross hairs show the most likely value in the two-dimensional projection of the full probability distribution function. 



nificance level. Furthermore, many of the best-fit single- 
component models favor stellar population ages that are older 
th an the age of the Universe for the measured redshift (see 
§ I6.1.2t . Restricting these models to ages less than the 
age of the Universe increases the minimum x" value, mak- 
ing the difference between single- and two-component fits 
more pronounced. We interpret this behavior to indicate that 
in general the DRG population has star-formation histories 
that are more complex than simple monotonically evolving 
stellar populations. A similar scenario has been suggested 
base d on m odeling the SEDs of LBGs at these redshifts 
(e.g.,ISawicki & Yee 19 98; Papovich. Dickinson. & Fergusoij 
I29OI; Panovic h et alJ l2004ar Shanlev et al. 20 0j), and is 
likely consiste nt with hier archical models dDe Lucia et alJ 
12005: Nagamine et alJ2005h . 

6. DISCUSSION 

The ensemble properties of the DRGs span a range of the 
stellar-population model parameter space. Broadly speaking, 
the (/-^L5)vega > 2.3 mag color selection identifies galaxies 
whose rest-frame optical and near-IR light is dominated by 



later-type stars, and galaxies whose light is dominated by 
heavily extincted starbursts (F ranx et alJl2003ft . In this way, 
it is not dissimilar to traditional R- K or I - K selection cri- 
teria for EROs, but it tends to pick out objects at higher red- 
shifts. In addition, roughly ~ 15% of the DRGs are luminous 
in X-rays, implying that some of the UV to IR emission may 
stem from SMBH accretion processes rather than star forma- 
tion. Here we consider the implications that our analysis has 
for star-formation in massive galaxies at high redshifts. We 
briefly consider AGN activity in these galaxies, with a more 
detailed analysis to be presented in L. A. Moustakas et al. (in 
preparation). 

6.1. Stellar Populations and Star-Formation Histories of 

DRGs 

6.1.1. Dust Extinction 

The best-fitting single-component DRG models span ex- 
tinctions of ^0-6 mag at 1600 A. The mean dust ex- 
tinction from the best fits to the ACS to IRAC photom- 
etry is (A1600) =3.1 mag, with a standard deviation of 
1.5 mag. The mean is substantially larger than that inferred 
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Fig. 10. — Compai'ison of the minimum reduced from fitting tlie DRGs 
with models of different star-formation histories. The abscissa shows the 
minimum reduced derived for models with single-component, exponen- 
tially decaying star-formation histories. The ordinate shows the minimum re- 
duced X' for the two-component models where one component corresponds 
to the exponentially decaying star-formation histories and the other corre- 
sponds to an instantaneous burst of star-formation at Zfomi = oo. that evolves 
passively thereafter 



for UV-selected LBGs at comparable redshifts^ for which 

tiDovich. Dickinson. & Ferguson (200ll) and IShaplev et al. 
001) find average extinction values of 1.2 and 1.6 mag 
at 1600 A, respectively. Forster-Schreiber et alJ i^Q4) find 
that DRGs in the FIRES fields have median extinction of 
2± 6.2 mag at 1600 A for the Calzetti extinction law (averaged 
by number), but restricted to models with constant SFR, and 
these values chang e little after including Spitzer/IRAC data 
jLabbe et alJl20()l . Although on the surface this appears at 
odds with our findings, Forster-Schreiber et al. found that us- 
ing exponentially declining models reduces the measured ex- 
tinction by roughly a factor of two, in better agreement with 
our results. 

To study differences between the sub-populations of DRGs 
(i.e., those with substantial mass in late-type stars and those 
dominated by highly dust-extincted starbursts), we divided 
the sample into subsets with E{B-V) < 0.35 and E(B-V) > 
0.35. We have used the results from the fits to the single- 
component models for this selection, but there is little change 
if we use the two-component models. This division point 
is the approximate upper bound on values inferred for UV- 
selected LBGs (largely as a result of the UV-selection itself; 
e.g., Adelberger & Steidel 2000). 

DRGs with E(B-V) > 0.35 comprise roughly 40% 
(61/152) of the full sample. These galaxies have a mean red- 
shift of (z) ~ 1 .7 with standard deviation a(z) = 0.5. The mean 
dust extinction is 4.6 mag at 1600 A for the Calzetti law. We 
estimate the resulting IR luminosities from these extinction 
values using the empirical relations bet ween the UV lumi- 
nosity, extinction, a nd far-lR luminosity (tMeurer et alJll999t 
ICalzetti et al] |2000l). In these galaxies, the observed UV lu- 
minosities and model extinction (using the 68% confidence 
limits) correspond to IR luminosities of lo""'^^ Lq. These 
are generally less than the 24 /im-derived IR luminosity or the 
24 //m upper limit for undetected galaxies). However, because 



the 68% confidence interval on galaxy extinction is large, 
(5A(1600) 1 - 2 mag, in roughly --50% of the objects the 
extinction from the upper 68% confidence bound yields a Lir 
comparable or exceeding the 24 fim value. This is consistent 
with the fraction of DRGs in figure |5l with IR/UV rati os and 
UV spectral slopes near the local relation from Meurer et alJ 
^999). 

The DRGs with best-fit E(,B-V) < 0.35 have a mean ex- 
tinction of 2.0 mag at 1600 A. Later-type stars dominate the 
rest-frame optical and near-lR light in these galaxies, while 
the UV rest-frame emission stems from small amounts (by 
mass) of ongoing star-formation with low to moderate extinc- 
tion. The mean redshift for the DRGs with E(B-V) < 0.35 is 
(z) ~ 2.5, noticeably larger than that for the DRGs with higher 
extinction. This higher redshift range arises because galaxies 
with a strong Balmer/4000 A break satisfy the (y-A'j)vega > 
2.3 mag selection only for z > 2.0 as this b reak is moving 
through the J and H bands JFranx et alJ2003l) . whereas heav- 
ily reddened galaxies can enter the sample at lower redshifts. 
The lower E{B-V) estimates from the model fitting would 
smaller IR luminosities relative to the DRGs with higher ex- 
tinction described above. Unlike the case for the higher- 
extinction DRGs, the IR-luminosity for the the DRGs with 
E(B-V) < 0.35 are less than the 24 /im-derived Lir by fac- 
tors of w 2 - 20. Therefore, the observed 24 /im emission in 
DRGs with relatively low dust extinction from the SED mod- 
eling does not originate only from the extincted stellar pop- 
ulations that dominate the UV and optical rest-frame light. 
These galaxies require either additional embedded star forma- 
tion, or an obscured AGN, or both, and these obscured com- 
ponents contribute negligibly at bluer wavelengths. 

6.1.2. Stellar Population Ages and Star Formation Histories 

Our stellar-population modeling of the DRGs is more sen- 
sitive to the ratio of model age, f, to the star-formation e- 
folding time, r, than on f or t individually, and we dis- 
cuss them simultaneously. From the models with a sin- 
gle, exponentially decaying SFR the median age is 1.1 Gyr, 
and on average a DRG has undergone f/r ~ 4 e-folding 
times. In an analysis of the 34 DRGs in the FIRES fields, 
Forster-Schreiber et al. (2004) tested constant star formation 
histories (r ^ t) and found median stellar population ages 
of f ~ 1.7-2 Gyr The lower median ages for our models re- 
sult because they favor star-formation histories with a smaller 
number of f /r e-folding times than the constant star forma- 
tion histories. Indeed, Forster-Schreiber et al. also considered 
models with r = 300 Myr, for which they derive a median age 
of 1 Gyr, implying f/r ^ 3. This is consistent with our result 
taken over a wider range of star-formation histories. 

There are 17 DRGs (~ 11% of the sample) where the 
lower limits from the 68% confidence range from the single- 
component-model fits are older than the Universe. Of these, 
nine (53%) show indications of AGN activity either because 
they are detected in the Chandra data, they have Lir > 
10'^ Lq, and/or they have have IRAC colors satisfying the 
AGN selection of Stern et al. (2005, see also § l6.3t . We sus- 
pect an AGN may influence the observed photometry and the 
model fits. The remaining DRGs have SEDs qualitatively 
similar to the example shown in Figure but the model fits 
favor older ages. In all these cases the two-component model 
fits have solutions with ages less than the age of the Universe 
within the 68% confidence interval, and are thus more phys- 
ical. Moreover the reduced minimum values for the two- 
component fits are significantly improved for more than half 
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of these, implying such models are better realizations of the 
data. 

Thus, including the constraint that the age of the model 
stellar population be less than that of the Universe favors 
models with more complicated star-formation histories than 
the single-component, monotonically evolving models. This 
reinforces our conclusion that the two-component models 
better describe the star-formation histories of the DRGs 
based on their lower reduced Xo values (see § 15. 3> . The 
real star-formation histories are probably even more com- 
plex. It is unlikely that the star-formation histories of DRGs 
are consistent with scenarios where all their mass formed 
in the distant past with subsequent passive evolution (see, 
e.g. ICimatti et al. 2002b). Massive galaxies at these red- 
shifts continue forming stars up to the epoch we are observ- 
ing them. Recent semi-analytic prescriptions of the star- 
formation histories of massive galaxies typically involve both 
stochastic "burst" modes, and quasi-continuous "quiescent" 
modes of star formati on (e. ^., Somerville, Primack, & Faber 
I200U iDe Lucia et al. 2005: Nagamine et al. 2005), qualita- 
tively consistent with our analysis. It seems that the earliest 
star-formation episodes in DRGs did occur in the distant past, 
and that much of their stellar mass was assembled much ear- 
lier than the epoch at which we are observing them. Taking 
the median age of the stellar population as a fiducial value 
for the onset of star formation, DRGs began forming stars at 
Zform ^3.5. Their progenitors may be the star-forming galax- 
ies observed at those earlier epochs as UV-dropouts at z > 3.5 
(e.g.,^apovich et al. 2004a). 

6.L3. Stellar Masses 

The SED modeling provides relatively robust estimates of 
the stellar-mass-to-light ratios of these galaxies, providing 
good measures of their stellar masses. Figure fTTI shows the 
distribution of the stellar masses for the single- and two- 
component star-formation histories for both the DRGs and 
the galaxies in the HDF-N. Typical errors are 0. 1-0.3 dex (for 
a given IMF and metallicity). 

The median stellar masses for the GOODS-S DRGs are 1 . 1 
and 1.7 x 10" A4q for the single- and two-component star- 
formation history models, respectively. There is little differ- 
ence between the stellar masses of DRGs detected and unde- 
tected at 24 jim. The interquartile range (containing the inner 
50% of the sample) spans from 0.13 to 3.3 x 10" Mq for the 
single-component models, and from 0.29 to 4.6 x 10" A4q 
for the two-component models. The main effect of the 
two-component models is to increase the masses of objects 
with low masses from the one-component models. Higher 
mass objects typically have larger stellar-mass-to-light ratios, 
A4/L, and are less effected because the uncertainties onA4/L 
are smaller (e.g., Dickinson et al. 2003). For the 34 DRGs in 
the FIRES fields, iForster-Schreiber et alJ li2004ft found me- 
dian stellar masses of = 0.8-1.6 x lO" Mq, with an in- 
terquartile range of ~ 0.6-3 x lO" A4q (using different as- 
sumptions for the star-formation history). Because our mod- 
els are taken over a wider range of star-formation histories, 
we expect our stellar masses to be slightly lower with respect 
to those derived for constant star-formation models (Forster- 
Schreiber et al. 2004). Thus the range and average of the DRG 
stellar-mass distribution seem generally consistent with those 
in the FIRES fields (see also Labbe et al. 2005). The tail of 
lower-mass DRGs in our sample probably also arises from 
the greater number of objects with lower redshifts. Restrict- 
ing our sample to a higher minimum redshift would further 



increase the lower-bound of our interquartile range. 

The long-dashed line in Figure ITTI shows the stellar mass 
limit for a passively evolving stellar population formed in a 
single burst at z = cxd and with K, = 23.2 mag. Our DRG 
selection is approximately complete in stellar mass for pas- 
sively evolving galaxies with A4 > 10" A4q to 2 < z < 3, 
because no galaxy can have a mass-to-light ratio higher 
than a maximally old stellar population. Bluer galaxies 
with Ai > 10" A4q and lower mass-to-light ratios such as 
UV-luminous LBGs could be excluded by the (/-/rv)vega > 
2.3 mag selection. However, Shapley et al. ( 2005J infer that 
the majority of LBGs at z 2 with M > 10" A^f7^ woul d 
satisfy (/-/Ovega > 2.3 mag (see also. iReddv et al.ll2005ft . 
At z < 2, the J-Ks color selection may miss some passively 
evolving galaxies as the Balmer/4000 A break shifts to wave- 
length s below the blue-edge of the 7-band filter ( Daddi et alJ 
12004ft . We consider the DRG selection roughly complete in 
stellar mass to > 10" A4q over the majority of the red- 
shift range 2 < z < 3. 

Based on the fits to the single-component models, the inte- 
grated DRG stellar mass density for objects at 2 < z < 3 with 
M > 10" Mq is 2.4-^1 X 10^ Mq Mpc'^ This increases 
to 2.8;!;o5 ^ -^G) Mpc"-' using the stellar masses derived 
from the two-component model fits. In both cases the un- 
certainties are estimated using a bootstrap resampling of the 
dataset, which constructs random samples of DRGs with the 
sample size taken from the Poisson, counting uncertainties, 
and stellar masses drawn from the measured distribution (with 
replacement) modulated by the inferred stellar-mass errors. 
However, many of the DRGs show evidence for an AGN, ei- 
ther based on X-ray detections , IR luminosity, or rest-frame 
near- and mid-IR colors (see § l6.3t . These objects have rela- 
tively higher derived stellar masses (Figure [TTl. which could 
imply that AGN tend to reside in the most massive galaxies 
at these redshifts. Alternatively it could be that AGN may 
contribute to the rest-frame UV-to-near-IR emission, lead- 
ing to unduly larger stellar mass estimates. To bound this ef- 
fect, we recalculate the stellar mass densities excluding all 
DRGs with X-ray detections, Lir > 10'^ Lq, or with AGN- 
like rest-frame near-IR colors. In this case, the stellar mass 
density decreases to 2.0:l:o4 ^^'^ '^■^^'t ^ -^Q Mpc"^ for 
the single- and two-component models, respectively. In this 
case, the DRGs contribute 25-70% to the total stellar mass 
density integrated over all galaxies at 2 < z < 3, in reasonable 
agreement with the conclusions of iRudnick et alJ J2003I) and 
iFontana etaLc2003,) . 

6.1.4. Star-Formation Rates and "Dead" Objects 

The 24 ^m-detected DRGs span SFRs from - 100- 
1000 Mq yr"', excluding those objects directly detected in 
the Chandra data (see Figure |6jl. The mean SFR for these 
sources is ~ 500 Mq yr"', with a systematic uncertainty in 
the 24 /im to SFR conversion of « 0.5 dex. This value in- 
cludes objects with Lir > 10'^ Mq or rest-frame near-IR 
colors indicative of AGN (see § 16. 3> . Excluding these ob- 
jects, the average SFR of the 24 /^m-detected DRGs drops 
to ~ 400 Al© yr"'. The mean SFR for the complete DRG 
sample (including those not detected with MIPS at 24 ^m) 
is lower still. Taking the conservative limit that the 24 ^m- 
undetected DRGs have no ongoing star formation, the mean 
SFR is ~ 220 Mq yr"' excluding the X-ray sources. This 
does not change if we also exclude those sources with infrared 
luminosities or colors indicative of AGN. 
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Fig. 1 1 . — Stellar masses derived from modeling the galaxy SEDs, plotted as a function of redshift. The left panel shows the masses derived using the single- 
component, monotonically evolving star-formation histories. The right panel shows the masses derived using two-component models. In both panels filled red 
circles denote the GOODS-S DRGs undetected at 24 fj,m. and open red circles show those DRGs with 24 /im detections. Blue squares denote the galaxies in 
the HDF-N. Black stars show those DRGs with X-ray detections. Error bars are not shown for clarity, but mean error bars are shown as a function of galaxy 
stellar mass. The short-dashed line shows the value of the characteristic stellar mass of the present-day mass function (.Cole et al. .2Q01) . In the right panel, the 
long-dashed line shows the stellar-mass for a passively evolving stellar population formed instantaneously at z = 00 with Ks = 23.2, which is the stellar-mass 
limit for the GOODS-S A'j-band data for such systems. 



The mean SFR for the DRGs reported here is some- 
what larger, but comparable to recent measurements of the 
"stacked" X -ray a nd sub-mm emission from DRG samples. 
iRubin et alJ J2004I) and .Reddv etalJ ^05f) find the average 
SFR is ~ 100-300 Mq yr"' based on statistical X-ray detec- 
tions for DRGs with Ks < 23 AB mag. The inferred mean SFR 
for the various DRG populations is strongly dependent on the 
limiting /iTj-band of the survey, and in the case of the X-ray 
derived measurements, on indications of AGN activity within 
the sample (see the discussion in Reddy et al. 2005). How- 
ever, the X-ray SFR calibration has a systematic unce rtainty 
on the order of a factor of five (Ranalli e t al. 2003; Per sic et al.l 
120041) . owing to assumptions on the formation timescales of 
X-ray binaries in starbursts. A higher SFR-to-X-ray lumi- 
nosity calibration may be appropriate for galaxies in the more 
intense starbursts such as those for the DRGs ( Pers ic et al.l 
|2P04; Tens et al. 2005), so some scatt e r is ex pected in the 
conversion is expected. iKnudsen et alJ i2005l) report an av- 
erage SFR of 130 Mq yr~\ based on a stacked sub-mm 
850 fim flux density of DRGs with K, < 24.4 AB mag, but 
is sensitive to the assumed average dust temperature (e.g., 
[chapman et al. 2005). Their sample extends to DRGs roughly 
a magnitude deeper in the /Tj-band, and Reddv et al. (2005) 
demonstrate that the mean SFR of all types of star-forming 
galaxies at z ~ 2 (including DRGs) decreases with decreas- 
ing /Tj-band flux density (see also Daddi et al. 2005b). Given 
the systematic uncertainties in the SFR calibrations, and the 
varying limiting magnitudes of the different DRG samples, 
we conclude that the mean SFR we derive using the MIPS 
24 /im data for the GOODS-S DRGs is in broad agreement 
with these other values. The "typical" DRG is forming stars 
at rates in excess of > 200 A4q yr"'. 

We can set an upper limit on the lifetime of a starburst in 
the DRGs at the gas-consumption timescale, defined to be 



the ratio of the gas mass to the SFR. For a massive galaxy 
with a molecular gas mass of A4{H2) < 10" TW© (compara- 
ble to the gas reservoirs of the massive z > 3 radio galaxies 
'De Br euck e'tani2003l l2005t iGreve. Ivison. & Papadopoulosl 
2004), the gas-consumption timescale is < 10^- 10^ yr for 
SFRs of ~ 10^-10-^ Mq yr"' (although these timescales 
would be shorter for lower gas masses, such as those in lo- 
cal ULlRGs, see Downes & Solomon 1998). This timescale 
is generally shorter than or comparable to the median ages we 
derive for the DRGs (consistent also with findings in Forster- 
Schreiber et al. 2004). Recent theoretical work suggests that 
high-redshift, massive galaxies quasi-continuously form stars 
at high rates for periods of 1-2 Gyr, with shorter periods 
(< 100 Myr) of boosted star formation (Nagamine et al. 2005; 
Finlator et al. 2005; De Lucia et al. 2005). If the DRGs sus- 
tain these the high SFRs for more than ~ 10'^- 10^ yr, then 
they will exhaust their gas supply unless fresh cold gas is re- 
peatedly or continuously accreted. Keres et al. (2005) predict 
that cold-gas accretion may dominate in massive galaxies at 
high redshifts, which may advocate such a scenario exists for 
the massive DRGs. 

Assuming a lifetime of 10**- 10'^ yr, a galaxy with a SFR 
of 1O^-1O^A^0 yr"' would assemble 10'°"'^ Mq in stars. 
If LBGs are forming stars rapidly and/or experiencing recur- 
rent bursts at redshifts higher than that of the DRGs, z > 3.5, 
then eventually they will assemble a sufficient population of 
late-type stars to dominate the optical to near-lR rest-frame 
light. It thus takes of order 1 Gyr to produce a DRG (un- 
less high-redshift starbursts have an IMF weighted towards 
higher-mass stars, in which case even more time is required; 
Ferguson et al. 2002). Under these star-formation histories, 
the first progenitors of DRGs at z ~ 2 to 3 should appear as 
starburst galaxies at z > 3 to 5. Massive star-forming galaxies 
appear to exist by z ~ 5 from U V-dropout surveys ( Evles et alJ 
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than the 50% completeness limit, f,y(24fim) < 60 /iJy. We 
also restrict the DRG and HDF-N samples to those objects 
with 1.5 < z < 3.0, where we are approximately complete 
for massive galaxies. Furthermore, restricting the sample to 
z < 3 removes few galaxies from the samples, an d faci litates 
the contrast with the lower-redshift samples (see § I6.4> . 

The DRGs have higher stellar masses and SFRs relative to 
the HDF-N 24 /im-selected sample, owing primarily to the 
fainter flux limit of the 24 //m data in the HDF-N. There is 
a weak correlation between SFR and stellar mass, although 
this trend may be softened if the SFRs of DRGs with upper 
limits are in fact lower Objects with stellar masses greater 
than 10" J^Q span 1.5 dex in SFR, excluding the X-ray- 
detected DRGs with inferred Lir > 10'^ Lq. The upper en- 
velope of SFRs decreases for lower mass galaxies, and such 
objects should be detected in our data if present. Therefore, 
at these redshifts galaxies with the highest stellar masses also 
appear to be experiencing the highest SFRs. 



Fig. 12. — The SFRs derived from the .sum of the UV and IR luminosities 
plotted as a function of galaxy stellar mass. Red circles denote the GOODS- 
S DRGs with 24 fim detections; gray triangles show upper limits for DRGs 
undetected at 24 /im, assuming they have f,^{24fim) < 60 /ijy. Black stars 
indicate those DRGs with X-ray detections. Blue squares con'espond to the 
24 /^m-selected galaxies in the HDF-N. All galaxy samples in this plot are 
restricted to the redshift range 1.5 < z < 3.0. Systematics dominate the un- 
certainties on the SFRs and are fa 0.5 dex, as indicated by the inset error bars. 
The error bars also show the mean uncertainty on the stellar masses in three 
mass bins. 



I2005tf^n et alJ200^ . and such objects could be forming the 
stellar populations that evolve to become the old components 
at the redshifts of the DRGs. 

Based on the model fits, we find several candidates among 
the DRGs for passively evolving galaxies at high redshifts 
(see also §|3|l. There are 15 DRGs with best-fit models that 
favor older ages (t > 1 Gyr) that have undergone several e- 
folding times (f/r > 3) with low dust extinction {E[B-V] < 
0.1) and no X-ray or 24 /im emission. Massive, "dead" 
DRGs are either galaxies that are between starbursts but al- 
ready with a massive old stellar population, or galaxies that 
have completed the intense phase of their assembly. These 
galaxies span 1.9 < z < 3.0, and have ijj^ -Ks > 2.6 and 
Kg- [4.5/im] < 1.2 (see also Figure|31and discussion in §1^. 
These objects contribute ~4.9±1.4 x 10^ Mq Mpc"^ to 
the global stellar mass density at these redshifts. Such mas- 
sive, "dead" objects make up only a fairly low fraction of the 
stellar mass density (< 10%) integrated over all galaxies at 
z 2-3, and therefore are relatively rare at these redshifts. 
In contrast, early-type, massive (> lO" Mq) galaxies ac- 
count for ^ 30% of the total local stellar-mass density for the 
stellar-mass function of Bell et al. (2003, see also Baldry et 
al. 2004). Thus, the fraction of the global stellar-mass density 
in such galaxies increases threefold from z > 2 to 0. 

6.2. The Relation between Star Formation and Stellar Mass 
in Galaxies at Z'^ 1.5 — 3.0 

Figure El shows the SFRs as a function of the stellar 
mass for the DRGs and the HDF-N sample. The figure 
shows the stellar masses derived from the models with single- 
component star-formation histories, although using the two- 
component fits does not strongly change any of the conclu- 
sions. DRGs not detected at 24 /.tm by MIPS are shown as 
upper limits assuming that their 24 /im flux density is less 



6.3. Growth of Supermassive Black Holes in Massive, 
High-Redshift Galaxies 

Many of the DRGs with the highest stellar masses and 
SFRs are detected in the deep X-ray data. Twelve of the 
13 X-ray detected DRGs with M* > 10" Mq have X- 
ray to optical flux ratios, log/x//R > -1.0 (the remaining 
one has log/x//R ~ -1-4), and these imply the presence of 
an AGN with Lx > 10"*^ e rg s"^ for DRGs at z > 1.5 (e.g., 
iHornschemeier et al. 200 ll lAlexanderet 111 12003). In addi- 
tion, Worslev et al. (2005) conclude that as much as 40% 
of the hard X-ray background stems from heavily obscured 
AGN at z > 0.5. Recent Spitzer IR observations are find- 
ing many AGN candidates whose X-ray, UV, and optical 
emission is heavily obscured by gas and dust, and missed 
in deep X-ray surveys (e.g., Alonso-Herrero et al. '200^ 
Donlevetal. 2005). Alonso-Herrero et al. (2005) find that 
as many as 50% of MIPS 24 fj,m sources with red IRAC 
colors are not detected in the deep X-ray data. This sug- 
gests that AGN may be hidden behind sufficiently Compton- 
thick material that none of the direct X-ray emission escapes 
(Brandt & Hasinger 2005). Models of the X-ray background 
support this possibility (Gilli 2004). 

At 1.5 < z < 3.5, galaxies dominated by the light from 
stellar photospheres should have relatively blue [4.5/im] - 
[8.0/im] colors (see §|3}. Alternatively galaxies dominated by 
AGN emission in the rest-frame near-IR should have redder 
IRAC colors (see , ej., Rieke 1978; Neugebauer et al. 1979t 
i Lacv et all 12004 lAlonso-Herrero et al.. 2QQ^ Stern et all 
I2005h . iLacvetal J 120041) found that^vfert 1 galaxies and 
obscured AGN have red [4.5/im] - [8.0/im] colors, similar to 
the red rest-frame J-K, colors of Parkes Quasars in 2MASS 
at z ^ 0-0.5 (Francis et al. 2004). Similarly, of the non- 
X-ray detected DRGs with IRAC counterparts, roughl y 25% 
(26/106) satisfy the IR color selection for AGN of Ster n et alJ 
Ii2005ft . Upon inspection of their individual SEDs, most 
are plausible AGN candidates with red IRAC colors. Four 
of these DRGs are star-forming objects at z > 2.9 whose 
IRAC colors mimic those of AGN candidates at lower red- 
shift. These galaxies are optically fainter and have higher 
redshifts than typical galaxies in the Stern et al. sample, and 
are an additional source of contamination in IR selection of 
AGN. An additional 10 of these DRGs have z ~ 1.0-1.4 
and their SEDs are consistent with heavily extincted starbursts 
whose [8.0/im] magnitudes may be augmented by the 3 fim 
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PAH emission feature.'^ Excluding these 14 objects, approx- 
imately 10% (10/125) of the X-ray-undetected DRGs may 
have dust-enshrouded AGN that contribute some fraction of 
the rest-frame near- and mid-IR light. Combined with the 
15% of X-ray-detected DRGs, perhaps a quarter of the DRG 
population host AGN, which implies they are actively grow- 
ing their SMBHs. 

The X-ray-detected fraction of DRGs apparently rises with 
increasing IR luminosity (>50% at Lir > lO'-' Lq), a trend 
also seen in the analysis of X-ray emission from high-redshift 
sub-mm galaxies (Alexander et al. 2005). The GOODS-S 
field lacks deep sub-mm observations, which would allow us 
to explore any direct overlap between these coeval popula- 
tions. Locally, the most luminous ULIRGs and HyLlRGs fre- 
quently host AGN sufficiently powerful to affect their global 
properties ( Veilleux et al. 1995). If the IR-luminous DRGs 
are similar, then the high incidence of AGN is expected, es- 
pecially for those with inferred Lir > 10'^ Lq. This implies 
that the massive galaxies at these epochs may be both assem- 
bling their stellar populations and growing SMBHs (see also 
L. Moustakas et al. in preparation). Locally, the space density 
of ULIRGs is w lO'^ho Mpc"^ (Sanders et al. 2003), while 
the DRGs with ULIRG-like luminosities have a space density 
of 5.8 ±0.9 X 10~^/i7o Mpc"^ — an increase by a factor of 
600 (and similar to the conclusion for Bz/T-selected objects, 
^addi et al. 2005b). Thus, such objects are a much more com- 
mon phenomenon at high redshift. 

6.4. Evolution of the Star Formation Rate as Function of 
Stellar Mass 

The DRGs have higher specific SFRs (SFR per unit stellar 
mass, than galaxies with comparable mass at lower 

redshifts (z < 1). Figure [T3l shows the specific SFRs for 
DRGs and HDF-N galaxies at 1 .5 < z < 3.0, comparing them 
t o gal axies at lower redshift in the COMBO-17 survey (see 
The SFRs for the COMBO-17 galaxies are calculated 
using the MIPS 24 /im imaging and rest-frame UV emis- 
sion in an analogous manner as for the DRGs. Masses for 
COMBO-17 galaxies were estimated from their rest-frame 
M(V) and U -V colors, and have a typical uncertainty of 
0.3 dex. This method will tend to overestimate the masses of 
galaxies involved in starbursts, but it is fairly robust for qui- 
escent galaxies such as those that dominate our conclusions 
here (see Bell et al. 2004, 2005a). 

We are biased against galaxies with low stellar masses and 
low specific SFRs, which causes the lower "envelope" in all 
the panels of Figure [O] This selection effect can produce 
the apparent anticorrelation between specific SFR and stellar 
mass (most apparent in the COMBO-17 galaxy plots). The 
galaxies at z < 1 show an upper envelope in the sense that 
there is a lack of galaxies with high specific SFRs and high 
stellar masses. In contrast, the DRGs at 1.5 < z < 3 show a 
nearly unchanging range of specific SFRs at all stellar masses, 
although there is a hint that the DRGs at the high-mass end 
(A4 > lO" Aio) show an upper envelope on their specific 
SFRs. Even then the massive galaxies have systematically 
higher specific SFRs than the z < 1 galaxies in COMBO-17. 
Quantitatively, the DRGs with X > 10" A^q and 1 .5 < z < 3 
have specific SFRs of 'i' / A4 ^ 0.2- 10 Gyr"', with a mean 

Stern et al. (2005) find that 17% of objects satisfying their IRAC color- 
color selection are spectroscopically classified as galaxies, roughly consistent 
with the galaxy contamination here, modulo differences in survey limiting 
magnitudes. 



value of ^ 2.4 Gyr ' (excluding those with X-ray detections). 
By Z--0.7 galaxies with X > 10" Mq have '^/M^O.l- 
1 Gyr"' and at z 0.4 galaxies with A< > 10" Mq have 
^t/A^ < 0.5 Gyr"' — an order of magnitude lower than for 
the massive DRGs. 

This downward evolution in the specific SFRs of massive 
galaxies has been referred to as "downsizing" (see §0. The 
hypothesis is that massive galaxies host most of the SFR den- 
sity at high redshifts, and that galaxy formation shifts to less- 
massive systems at lower redshifts. It is unclear if downsizing 
is a proper description as nearly all galaxies at 1.5 < z < 3.0 
have higher specific SFRs at higher redshift, not just the most 
massive. "Downsizing" is not equivalent to a global change in 
'i'/Ai with redshift. However, our results indicate that star- 
formation in massive galaxies is reduced for z < 1 as galaxies 
with lower stellar masses have higher specific SFRs. It ap- 
pears that massive galaxies are largely done assembling their 
stellar mass in high specific-SFR events by z ~ 1 .5 (although 
mergers accompanied by low specific SFRs, i.e., "dry" merg- 
ers, could still take place at lower redshift, see, e.g.. Bell et al. 
2005b; Faber et al. 2005; van Dokkum 2005). 

We define the integrated specific SFR, T, as the ratio of the 
sum of the SFRs, ^t,, to the sum of their stellar masses, Aii, 
over all / galaxies. 



This is essentially just the ratio of the SFR density to the 
stellar mass density for a volume-limited sample of galax- 
ies. Figure O shows the integrated specific SFRs for DRGs 
at 1.5 < z < 3.0 and COMBO-17 galaxies at z 0.4 and 

0. 7 with M > 10" Mq. The data point for the DRGs in- 
cludes all DRGs with M > 10" Mq, and assumes that the 
DRGs without 24 fim detections have a no star-formation 
(vp = Mq yr"'). The error box illustrates how changes 
in our assumptions affect the result. The bottom bound of 
the box shows the value if we exclude those objects with X- 
ray detectio ns, o r with IR luminosities or colors indicative of 
AGN (see § I6.3> . The upper bound of the box shows the value 
for DRGs in this redshift and mass range, and if we calculate 
SFRs for the DRGs without 24 /im detections assuming they 
have /^(24/tm) = 60 /iJy. Note that for this calculation we 
have used the stellar masses from the single-component fits. 
Using those from the two-component fits increases the DRG 
point and error box by a small amount, <0.05 dex. 

The integrated specific SFR in massive galaxies declines 
by more than an order of magnitude from z 1 .5 - 3 to red- 
shifts z ~ 0.7 and 0.4. Using the rest-frame UV luminosi- 
ties and [O II] emission -line fluxes o f galax ies in the Gem- 
ini Deep-Deep Survey, 'Junea u et al.l (l2005h found that the 
specific SFRs of galaxies with M > 10'"^ Mq decline by 
roughly a factor of six from z ~ 2 to 1 . Our results suggest 
that this value may be underestimated when the IR emission 
from these galaxies is included in the SFR determination. 

Figure [H] shows the specific SFR integrated over all 
galaxies by taking the ratio of the cosmic SFR den - 
sity to the integrated SFR density from ICole et ak (2001"), 

1. e., T = p*/ f p^dt and also from the predictions from 
Hernauist & Springel (2003). The integrated specific SFR 
over all galaxies declines steadily with decreasing redshift, 
because the SFR density peaks between z ~ 1 - 3 (Hopkii^ 
120041) . corresponding to when the stellar mass density appears 
to have grown the most rapidly CDickinson et al. 2003i) . That 
is, there is a decrease in the global specific SFR. 
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Fig. 13. — The specific SFR as a function of gala xy s tellar mass. The top panel shows the results for the DRG sample and the HDF-N galaxies, all restricted to 
1 .5 < z < 3.0, with symbols and errors as in Figure [l2l The bottom panels show the specific SFRs as a function of stellar mass for lower-redshift galaxy samples 
from the COMBO-17 survey (as labeled). Open circles show COMBO-17 galaxies detected with Spitzer at 24 fim, and small filled circles show upper limits for 
COMBO-17 galaxies undetected at 24 /xm. 



The evolution in the integrated specific SFR in massive 
galaxies is accelerated relative to the integrated global value 
for all galaxies. Galaxies with M > lO" Mq at z ^ 1.5-3 
were forming stars at or slightly above the rate integrated 
over all galaxies. These massive galaxies contribute between 
10-30% to the global SFR density at z 1.5-3 (compar- 
ing to Hopkins .2004.) . In contrast, by z < 1, galaxies with 
AA > lO'^ Mq have an integrated specific SFR much lower 
than value integrated over all galaxies. Thus, our results in- 
dicate that by z < 1 massive galaxies have formed most of 
their stellar mass, and lower-mass galaxies dominate the cos- 
mic SFR density. If we tried to track individual galaxies we 
would need to account for their growth in mass from z ^ 3 to 
z < 1 , and as a result this trend would be even stronger 

Our conclusion would be undermined if there were a sub- 
stantial population of passively evolving, massive galaxies 
at high redshifts that we have missed in our selection. The 
GOODS-S DRG sample is approximately mass-limited for 
passively evolving galaxies to 10" Mq at 2 < z < 3, al- 



though we are presumably missing some at z < 2 (see Daddi 
et al. 2004, and § I6.1.3> . However, if we restrict the calcu- 
lation to those DRGs with 2 < z < 3, the integrated specific 
SFR increases by roughly 0. 1 dex. Including those sources 
at z < 2 provides a more conservative measure of the evo- 
lution in the specific SFRs of massive galaxies, so we keep 
these objects here. The DRG selection is potentially bi- 
ased against other star-forming galaxies, most notably UV- 
luminous LBGs. However, LBGs with these stellar masses 
have high inferred SFRs, comparable even to those of the 
DRGs (see, e.g., Shapley et al. 2003; and Figure |5|i- More- 
over, Shapley et al. find that the inferred J-Ks colors of most 
UV-selected LBGs at z 2 with M > lO" Mq satisfy the 
DRG selection criterion. Therefore, including massive LBGs 
into the integrated specific SFR would have little effect on our 
conclusions, and we are confident that the results presented 
here apply to the general galaxy population at z ^ 1.5-3 with 
M > 10" Mq. 
This decline in the specific SFR of massive galaxies seems 
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Fig. 14. — Evolution of the integrated specific SFR, i.e., the ratio of the total SFR to the total stellar mass. The curves show the expected evolution of the ratio 
of the total SFR to the total galaxy stellar mass densities from an empirical fit to the evolution of the SFR density (solid lines, thick line includes correction for 
dust extinction; Cole et al. 2001), and the model of Hernauist & Sprinael (dashed line; 2003). The data points show results for galaxies with A4 > lO" A4q. 
The filled diamonds show the mean values derived for COMBO-17 galaxies, and the filled circle shows the mean value for the DRGs. The data point for the 
DRGs sets the SFR of those galaxies without MIPS detections to zero. The lower bound of the box around the DRG point shows the value if we exclude those 
DRGs with direct X-ray detections, and those DRGs with Ljr > lO'^^ Lq or IR colors indicative of AGN (see text). The upper bound shows the value for the 
DRGs with no MIPS detection have fu(2A^m) = 60 /iJy. The error bars on the data points themselves are derived by a bootstrap resampling of the dataset. These 
do not include systematic uncertainties in the SFRs, which are indicated by the inset error bar. 



to signify the point where massive galaxies have formed the 
bulk of their stellar populations. This may happen because 
they have converted most of their gas reservoirs into stellar 
mass, and thus can no longer sustain high SFRs. This may be 
supported by the timescales f or the galaxies to sustain their 
current SFRs discussed in § 16.1.41 Alternatively, because 
some of the massive DRGs show indications that they are fu- 
eling powerful AGN, it may be that AGN feedback provides 
a means to suppress further star formation in these objects 
((Hopkins et al. 2005). To understand the rapid evolution in 
the specific SFRs of massive galaxies at these redshifts, we 
will need to better understand the coeval assembly of stel- 
lar mass and AGN, and the energetic feedback from the pro- 
cesses. 

7. CONCLUSIONS 

We have investigated the properties of massive, star- 
forming galaxies at z ~ 1 - 3.5 using observations from the 



Hubble Space Telescope, ground-based near-infrared (IR) 
imaging, and IR observations from the Spitzer Space Tele- 
scope at 3-24 /im. From a /T^-selected galaxy sample over a 
~ 130 ai-cmin^ field in GOODS-S, we identified 153 DRGs 
with (7-/rs)vega > 2.3 to Ks < 23.2 mag, which is approxi- 
mately complete in stellar mass to lO" Mq for z < 3. The 
majority of DRGs (90%) are detected in deep Spitzer/IRAC 
imaging at 3.6-8.0 /im, and the remainder are confused 
with nearby sources within the IRAC beam. Roughly half 
of the DRGs are detected by SpitzerMlPS at 24 pm with 
fu{2A-pm) > 50 /iJy. At z 1-3.5, these 24 pm flux den- 
sities correspond to IR luminosities of Lir ^ 10"~'^ Lq. 

Based on the full suite of photometry, we find that the 
DRG selection identifies galaxies whose optical rest-frame 
light is dominated by a population of evolved stars combined 
with trace amounts of ongoing star formation (Zmed ^ 2.5), 
and galaxies whose light is dominated by heavily extincted 
(■A 1600 ^4-6 mag) starbursts (Zmed ^ 1-7), in roughly equal 
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proportions. Only a very small fraction of the DRGs (< 10%) 
in the sample have SEDs consistent with pure old (> 1 Gyr) 
stellar populations with no indication of current star forma- 
tion. We estimate that the number density of massive, pas- 
sively evolving DRGs at 2 < z < 3 is 3.2 x 10"^ /i^j] Mpc"\ 
which is nearly an order of magnitude lower than that for 
DRGs in the HDF-S (Labbe et al. 2005), although some of 
this discrepancy likely results from the fact that the HDF-S 
sample extends to fainter Ks-hand sources. The candidates for 
passively-evolving DRGs at 2 < z < 3 with M > 10" Mq 
contribute only a small fraction (< 10%) to the global stellar- 
mass density integrated over all galaxies at these redshifts. 
Thus, massive, passively evolving objects are rare at z > 2. 
However, the fraction of the stellar-mass density in massive, 
passively evolving galaxies increases threefold from z > 2 to 
0. 

The UV-derived SFRs for the DRGs are lower than those 
that include the reradiated IR emission by up to two orders 
of magnitude. The DRGs have IR/UV luminosity ratios typi- 
cally in excess of what is expected from their rest-frame UV 
spectral slope, consistent with observations of local LIRGs 
and ULIRGs. 

We compare the DRG photometry to stellar population syn- 
thesis models to estimate the stellar masses and to study 
the properties of the stars that dominate the rest-frame UV 
through near-IR light. Models allowing for a previous stel- 
lar population formed in an instantaneous burst in the distant 
past accompanied by additional on-going star-formation gen- 
erally provide better fits to the data. We conclude that the 
DRGs have complex and stochastic star-formation histories 
consistent with other star-forming galaxies at these redshifts 
and locally, and with predictions from hierarchical models. 

DRGs at z ~ 1.5-3 with stellar masses greater than 
lO" A4q have specific SFRs ranging from 0.2- 10 Gyr"'. 
This is more than an order of magnitude larger than that 
derived for galaxies with stellar masses above 10" Mq at 
z ~ 0.7 and 0.4. Simultaneously, the most luminous and mas- 
sive DRGs show indications for the presence of AGN either 
based on X-ray luminosity, IR luminosities, or IR colors. We 
find that as many as one-quarter of the DRG population con- 
tain AGN, and therefore the growth of SMBHs coincides with 
the formation of massive galaxies at z > 1 .5. Further implica- 



tions from the prevalence of AGN in DRGs will be discussed 
in L. Moustakas et al. (in preparation). 

At z ^ 1.5, galaxies with stellar masses > 10" AAq are 
forming stars at rates at or slightly higher than the global value 
integrated over all galaxies at this epoch. In contrast, galaxies 
at z < 0.7 with Ai > 10" Mq are forming stars at rates less 
than the global value for all galaxies. The evolution in the 
specific SFRs of massive galaxies occurs at an accelerated rate 
compared to that of all galaxies. The bulk of star formation 
in massive galaxies occurs at early cosmic epochs and that 
further mass assembly in these galaxies is accompanied by 
low specific SFRs. At z < 1 massive galaxies have formed 
most of their stellar mass, and lower mass galaxies dominate 
the cosmic SFR density. 
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